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Abstract 
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well  structures  were  studied  with  time-resolved  photoluminescence  utilizing  upconversion,  a  non-linear  wave 
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radiatively.  Einally,  spectral  data  was  taken  to  examine  the  spectral  decay  of  the  luminescence.  Carrier 
temperatures  were  extracted  from  the  spectral  data.  Type  I  structures  were  found  to  have  hotter  carrier 
temperatures  and  higher  Auger  coefficients  than  type  II  structures. 
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Abstract 

The  objective  of  this  researeh  is  to  improve  the  performance  of  semiconductor 
quantum  well  lasers  operating  in  the  mid-infrared.  Researeh  in  mid-infrared  laser  technology 
has  uncovered  numerous  applications  for  commercial  and  government  use.  Many  military 
applieations  require  operation  near  room  temperature  with  eompact  power  supplies.  Limiting 
faetors  for  improving  the  efficiency  of  mid-infrared  semiconductor  lasers  at  elevated 
temperatures  inelude  non-radiative  recombination  and  carrier  relaxation.  Carrier  lifetimes 
are  typically  on  the  order  of  5  nsec,  therefore,  reeombination  meehanisms  and  carrier 
relaxation  times  are  difficult  to  measure  directly.  Time-resolved  photo  luminescence  (TRPL), 
utilizing  a  non-linear  wave  mixing  teehnique  known  as  upeonversion,  is  used  to  study  the 
time  evolution  of  earrier  densities  in  semieonduetor  lasers.  Carrier  reeombination  and 
relaxation  information  is  extraeted  from  this  data  to  further  the  understanding  of  how 
quantum  well  struetures  impact  semiconductor  laser  performance. 

InAsSb/InAlASb  type  I  andInAs/GalnSb  type  II  multiple  quantum  well  structures 

3 1 

were  studied.  Auger  reeombination  coefficients  were  determined  to  be  3.3  +/-  0.1  x  10' 
em  /s  for  type  I  structures  and  6.0  +/-  1.0  x  10'  em  /s  for  type  II  struetures  at  77  K.  The 
Auger  coefficient,  whieh  describes  a  non-radiative  process  that  ean  dominate  when  using 
wavelengths  in  the  mid-IR,  was  found  to  be  higher  in  type  I  than  type  II  struetures.  This 
trend  is  in  agreement  with  previous  research,  however,  there  has  been  a  wide  range  of 
reported  values  for  the  Auger  eoefficient  in  similar  material.  Therefore,  careful  data 
reduetion  and  error  analysis  is  included  to  bound  the  reeombination  coefficients  with  some 
degree  of  certainty.  Carrier  temperatures  were  extraeted  from  the  speetral  TRPL  data. 


xii 


Type  I  structures  were  found  to  have  hotter  carrier  temperatures  than  type  II  structures, 
which  is  consistent  with  the  Auger  coefficient  results. 

TRPL  data  was  also  used  to  investigate  carrier  relaxation,  or  the  time  required  for 
energetic  carrier  will  cool  to  the  bottom  of  a  quantum  well,  where  it  can  recombine 
radiatively.  The  carrier  relaxation  times  were  found  to  increase  with  excitation  power  and 
decrease  as  the  quantum  well  width  increased.  Carrier  relaxation  times  for  the  type  II 
structures  were  16,  8,  and  5  ps  for  InAs  wells  of  width  20,  23,  and  26,  respectively.  These 
relaxation  times  were  measured  at  300  mW  excitation  power  with  a  fixed  GalnSb  well  width 
of  24  Angstroms. 


CARRIER  DYNAMICS  IN  MID-INFRARED  QUANTUM  WEEE  EASERS  USING 
TIME-RESOEVED  PHOTOEUMINESCENCE 

1.  Introduction 

Mid-infrared  (IR)  semiconductor  lasers  are  well  suited  for  a  variety  of 
applications  including  pollution  monitoring,  optical  communications,  and  laser  surgical 
equipment.  The  Air  Force  is  particularly  interested  in  applying  this  technology  to  IR 
counter-measures  and  laser  radar.  IR  guided  missiles  are  abundant  and  proliferated 
worldwide,  and  IR  passive  detection  systems  are  being  incorporated  into  many  fourth 
generation  aircraft.  Antimony  based  semiconductors  are  a  good  choice  for  IR  laser 
devices  because  of  their  flexibility  to  be  grown  to  take  advantage  of  the  atmospheric 
transmission  window  from  3-5  pm.  While  the  size  of  such  devices  is  attractive  for 
aircraft  applications,  the  efficiency  has  been  marginal.  This  paper  will  focus  on  two 
processes,  carrier  recombination  and  carrier  relaxation,  that  have  a  significant  effect  on 
the  performance  mid-IR  semiconductor  devices. 

1,1.  Motivation 

There  are  over  500,000  IR  guided  missiles  proliferated  worldwide  that  pose  a 
significant  threat  to  allied  aircraft  (Sawyer,  2000).  To  illustrate  this  threat.  Figure  1.1 
shows  the  number  of  aircraft  that  have  been  lost  to  IR  missiles  in  several  conflicts  dating 
back  to  1950’s.  Man  Portable  Air  Defense  Systems  (MANPADS)  are  IR  weapons  of 
particular  concern  to  mission  planners  because  of  their  mobility.  Avoidance  has  been  the 
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primary  tactic  to  defeat  the  MANPAD  threat  in  reeent  conflicts.  In  Operations  ALLIED 
FORCE  and  ENDURING  FREEDOM,  eombat  aireraft  in  hostile  territory  were  flown 
above  15,000  ft  with  few  exeeptions  (Lambeth,  2001).  Airdrops  and  visual  target 
identification  become  diffieult  at  this  altitude,  and  operators  are  further  restricted  by 
weather  eonstraints.  Unfortunately,  avoidanee  may  not  be  a  plausible  tactic.  When 
operating  out  of  foreign  airfields,  aircraft  pass  through  the  heart  of  the  MANPAD 
engagement  envelopes  on  take-off  and  landing.  In  an  effort  to  counter  the  IR  threat,  the 
Air  Force  has  researched  a  number  of  laser  based  IR  countermeasures  (IRCM). 
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Figure  1.1.  Aircraft  losses  due  to  IR  missiles  in  several  conflicts.  During  Operation 
Allied  Force  in  Kosovo,  RF  missiles  downed  two  aireraft,  IR  missiles  struek  two  aireraft, 
and  AAA  damaged  two  aircraft  (Kent,  2000  and  Lambeth,  2001). 
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The  Sensors  Directorate  of  the  Air  Force  Research  Laboratory  established  the  Laser 
Infrared  Flyout  Experiment  (LIFE)  program  to  test  advanced  IRCM  protection  of  large 
aircraft  (Taylor,  2001).  Eigure  1.2  (a)  shows  one  of  the  prototypes  for  the  EIEE  program 
developed  by  Sanders  (Sanders,  2000).  A  large  aircraft  would  carry  up  to  four  laser 
IRCM  systems  to  provide  all-aspect  protection.  Eighter  aircraft,  however,  lack  the  space 
and  power  required  to  operate  these  countermeasures.  Eurther  research  is  being 
conducted  to  reduce  the  size  and  improve  the  efficiency  of  laser  countermeasures  to  field 
an  IRCM  system  for  fighter  aircraft.  Eigure  1.2  (b)  shows  a  semiconductor  laser 
designed  to  achieve  that  goal. 


(a)  (b) 


Eigure  1.2.  Easer  IR  countermeasure  systems,  (a)  Advanced  Tactical  IRCM  System 
developed  by  Sanders.  The  insert  shows  the  rotating  mount  that  houses  the  laser  source 
(Sanders,  2000).  (b)  Easer  designed  for  IRCM  with  reduced  size  for  fighter  applications 
(Goyal  et  al.,  2000). 
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1.2.  Background 

Efforts  to  grow  quantum  well  laser  struetures  designed  to  operate  in  the  mid-IR 
range  have  foeused  on  InAsSb  type  I  struetures  (Turner,  1995;  Lane  et  ai,  1997).  Lane 
has  shown  that  these  struetures  require  half  the  threshold  eurrent  required  for  laser 
operation  than  similar  heterostrueture  lasers  (Lane  et  al,  1997).  Reeently,  type  II 
struetures  have  proven  even  more  effieient  by  deereasing  nonradiative  reeombination 
proeesses  through  band-gap  engineering  (Jang  and  Llatte,  1998).  Carrier  relaxation  has 
been  studied  in  type  II  materials  as  a  function  of  excitation  power  (Bogges  et  al,  2000), 
but  the  effect  of  varying  the  well  width  has  been  limited  to  single  quantum  wells  (Davis 
et  al,  1994).  Laser  performance  continues  to  improve  as  the  impact  of  the  growth 
process  is  better  understood. 

Time  resolved  photo  luminescence  (TRPL);  first  described  by  Mahr,  is  an 
experimental  method  that  provides  the  temporal  resolution  necessary  to  study  the  rapid 
recombination  of  carriers  in  semiconductor  devices  (Mahr  and  Hirsch,  1974).  The  TRPL 
technique  is  used  to  describe  the  evolution  of  carrier  densities  in  an  excited 
semiconductor.  Extensive  work  has  been  done  at  the  University  of  Iowa  to  determine 
recombination  rates  in  specially  grown  mid-IR  semiconductors  designed  to  reduce 
nonradiative  recombination  (Jang  and  Llatte,  1998;  Jang  et  al,  1996).  Work  at  ALIT  has 
been  done  to  determine  recombination  rates  (Cooley,  1997;  McKay,  2000),  and  carrier 
densities  based  on  spectral  width  (Lranz,  1997).  This  experiment  draws  upon  techniques 
and  methods  of  analysis  from  previous  work  in  an  effort  to  further  the  understanding  of 
carrier  dynamics  in  mid-IR  semiconductors. 
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1,3  Approach 

The  TRPL  experiment  uses  frequency  upconversion  to  study  the  characteristics  of 
photoluminescence  from  mid-IR  semiconductors.  Upconversion  is  a  method  of  wave 
mixing  used  to  measure  the  luminescence,  which  is  described  in  detail  by  Shah  (Shah, 
1988).  Multiple  methods  have  been  developed  to  detect  luminescence  over  very  short 
time  periods,  including  direct  detection  with  streak  cameras,  photoconductive  response, 
and  pump-probe  spectroscopy.  Upconversion  is  considered  the  best  technique  for 
studying  carrier  recombination  because  it  provides  the  best  temporal  resolution  and  does 
not  require  any  modifications  to  the  sample  (Jang  and  Flatte,  1998).  The  TRPL 
experiment  can  be  used  to  study  the  temporal  and  spectral  aspects  of  the  luminescence, 
both  of  which  were  accomplished  for  this  paper. 

The  TRPL  experiment  describes  the  photon  intensity  over  time.  To  extract 
meaningful  information  about  the  semiconductor  itself,  a  relationship  between  carrier 
density  and  luminescence  intensity  must  be  determined.  The  initial  carrier  density 
created  in  the  semiconductor  by  an  excitation  pulse  is  calculated  for  multiple  excitation 
powers  (Sermage,  1986).  These  theoretical  results  are  then  matched  to  experimental 
luminescence  intensities  at  the  same  excitation  power.  A  function  is  then  fit  to  these  data 
points  to  describe  carrier  density  in  terms  of  luminescence  intensity.  At  this  point,  the 
recombination  rate  equation  can  be  fit  to  the  carrier  density  data,  and  the  recombination 

The  recombination  rate  equation  is: 

2  3 

A  SRH  n+  B  j-ad'H  +  C  Auger'H 

(1.1) 


coefficients  can  be  extracted. 

dn 

dt 
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where  n  is  the  earrier  density  and  the  Asrh,  Brad,  CAuger  eoeffieients  deseribe  the  influenee 
that  a  partieular  reeombination  meehanism  has  on  the  earrier  density.  The  goal  is  to 
define  these  eoeffieients  for  different  semieonduetors  to  determine  struetural  effeets  on 
device  efficiency. 

The  initial  rise  of  the  carrier  density  can  also  be  isolated  to  examine  the  carrier 
relaxation.  The  rise  in  the  luminescence  signal  is  a  result  of  energetic  carriers  cooling  to 
the  bottom  of  the  quantum  well.  The  excitation  energy  and  well  width  were  varied  to 
study  their  impact  on  carrier  relaxation  time. 

The  TRPL  experiment  can  provide  spectral  information  about  the  luminescence  at 
fixed  points  in  time.  By  recording  data  at  a  series  of  times,  the  spectral  evolution  of  the 
luminescence  can  be  studied.  The  high-energy  tail  of  each  spectrum  was  used  to  extract 
the  effective  carrier  temperature.  The  evolution  of  carrier  temperature  provides  further 
insight  as  to  how  quickly  carriers  give  up  energy  to  the  lattice. 

1,4  Overview 

Chapter  2  will  discuss  the  theory  of  carrier  dynamics  in  a  quantum  well  structure. 
Chapter  3  provides  details  of  the  TRPL  experiment  and  the  upconversion  theory.  The 
results  and  data  analysis  are  described  in  Chapter  4.  The  carrier  densities  as  a  function  of 
time  for  four  samples  will  be  presented.  The  spectral  evolution  of  the  luminescence  from 
two  samples  will  also  be  shown.  Finally,  Chapter  5  will  discuss  conclusions  drawn  from 
this  analysis  and  recommendations  for  future  work. 
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2,  Theory 


Quantum  well  (QW)  structures  have  proved  useful  in  increasing  the  efficiency  of 
semiconductor  lasers.  A  lower  threshold  current  is  required  to  create  a  population 
inversion  in  QW  structures  because  they  provide  carrier  confinement  at  the  band  edge.  In 
addition,  there  is  evidence  that  nonradiative  processes  are  reduced  in  QW  structures 
(Hausser,  1990).  The  benefits  of  a  single  QW  structure  are  outweighed  by  the  small  gain 
medium  the  well  creates.  A  series  of  heterojunctions  with  alternating  band-gap  energies 
form  a  multiple  quantum  well  (MQW)  structure  with  a  much  wider  active  region  than  a 
single  quantum  well.  This  chapter  will  discuss  the  theory  that  dictates  the  carrier 
dynamics  in  a  MQW,  to  include  carrier  distributions,  recombination  mechanisms,  and 
scattering.  It  is  assumed  that  the  reader  has  a  basic  understanding  of  solid-state  physics 
and  is  familiar  with  carriers  in  bulk  semiconductor  materials. 


2.1,  Density  of  States 

Epitaxial  growth  processes  are  used  to  create  thin  layers  of  semiconductors  with 
varying  band-gap  energies.  If  the  thickness  of  the  narrow  band-gap  semiconductor  is 
similar  to  or  smaller  than  the  electron  de  Broglie  wavelength,  quantum  effects  become 
important  and  change  the  properties  of  the  semiconductor  (Shah,  1996).  The  narrow 
band-gap  semiconductor  dictates  the  well  width,  which  becomes  the  confined  region  for 
electrons  in  one  dimension.  The  electron  energy  in  an  infinite  quantum  well  is  quantized 
and  dependent  on  the  well  width  (L)  by: 


8  a  — 


1  2  2 
h  a 


-me-L 


(2.1) 
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where  Sa  is  the  energy  above  the  bottom  up  the  eonduetion  band,  a  is  an  integer,  h 
is  Planek’s  eonstant,  and  mg  is  the  eleetron  effeetive  mass.  While  the  quantum  well 
struetures  in  this  experiment  are  finite,  the  theory  of  quantized  energy  states  still  applies. 
The  caleulation  of  these  states  is  eomplex  and  beyond  the  seope  of  this  work.  The  3-D 
density  of  states  funetion  normally  applied  to  bulk  material  is  no  longer  valid  in  a 
quantum  well  beeause  earriers  are  now  eonfined  in  one  dimension.  The  density  of  states 
at  a  given  energy  level  in  the  2-D  and  3-D  eases  is  shown  below  (MeKelvey, 
1993:134,151): 

4-71  -me 

g2D  =  - ; — 

h 


g3D(E)  =  4-71 


2-mf 


I  J 


(2.2) 


The  2-D  density  of  states  is  independent  of  energy,  whieh  is  a  result  of  the 
quantized  energy  bands.  This  density  of  states  funetion  for  a  quantum  well  is  depieted  in 
Figure  2.1  (Saleh  and  Teieh,  1991:571).  The  density  of  states  funetion  will  beeome 
perturbed  at  high  exeitation  levels  as  a  result  of  eoulomb  sereening.  This  perturbation 
will  eause  a  narrowing  of  the  energy  band-gap  of  a  semieonduetor,  and  radiative 
reeombination  will  take  plaee  slightly  below  the  band-gap  energy.  This  eharaeteristie  is 
ealled  band  tailing  and  ean  be  seen  in  photolumineseenee  eurves  diseussed  later  in  this 
ehapter  and  experimental  data. 
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Density  of  States 


Figure  2.1.  Density  of  States  in  a  quantum  well  and  bulk  material.  The  dashed  line 
represents  a  3-D  free  partiele,  and  the  solid  line  represents  a  3-D  particle  confined  to  a  1- 
D  quantum  well.  Ea  and  Eb  are  the  energy  gaps  that  make  up  the  barriers  and  wells.  Ei 
and  E2  represent  the  first  two  confined  energy  states  in  the  well  of  width  E. 


2,2,  Carrier  Distribution  Function 

The  density  of  states  and  the  distribution  of  carrier  energy  dictate  the  number  of  carriers 
in  a  semiconductor.  The  Fermi-Dirac  distribution  function  describes  how  carriers  are 
spread  about  a  single  Fermi  energy  in  an  equilibrium  situation.  In  this  experiment,  a  non- 
equilibrium  situation  is  formed  when  a  pulsed  laser  creates  electron-hole  pairs.  Two 
separate  Fermi  functions  are  used  to  describe  the  electron  and  hole  distribution  about 
quasi-Fermi  energy  levels.  The  non-equilibrium  distribution  functions  for  electrons  (fi) 
and  holes  (fp)  are; 


fn(E) 


1  -I-  exp 


(E-Efh)' 


k-T 


l-fp(E)  = 


1  -I-  exp 


(E-Efh)' 


k-T 


(2.3) 
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where  £&  is  the  eleetron  quasi-Fermi  level  in  the  eonduetion  band,  Efp  is  the  quasi-Fermi 
level  in  the  valence  band,  and  T  is  the  carrier  temperature  (Bhattacharya,  1997:1 14).  The 
probability  of  a  recombination  event  occurring  is  the  product  of  the  electron  and  hole 
distribution  functions.  Assuming  that  the  excitation  source  has  been  removed  and  the 
quasi-Fermi  levels  are  deep  in  the  bands,  the  carrier  distribution  functions  can  be 
approximated  by  (Saleh  and  Teich,  1991:584): 

-E 

f(E)  =  A-e^'^  (2.4) 

The  carrier  distribution  and  the  density  of  states  together  describe  the  number  of 
carriers  available  for  recombination.  This  can  be  seen  in  the  notional  photoluminescence 
(PE)  curve  shown  below  in  Eigure  2.2. 


Eigure  2.2.  Example  PE  curves.  Low  energy  side  is  governed  by  density  of  states 
function.  High-energy  tail  is  due  to  carrier  distribution  function.  The  band  tailing  effect 
is  shown  as  energy  states  available  below  the  band-gap  energy  Eg. 
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The  density  of  states  funetion  dictates  the  low  energy  portion  of  the  PL  curve. 
The  high-energy  side  of  the  curve  is  a  result  of  carrier  energy  distribution,  and  will  vary 
with  temperature.  As  the  carrier  temperature  decreases,  the  high-energy  portion  of  a  PL 
curve  will  have  an  increased  negative  slope  and  the  width  of  the  curve  will  decrease. 

2.3.  Recombination  Mechanisms 

Radiative  recombination  in  semiconductors  is  measured  easily  and  can  be 
described  by  luminescence  curves  like  the  one  depicted  in  Figure  2.2.  When  electron- 
hole  pairs  are  created  with  an  excitation  pulse,  many  pairs  will  recombine  after  a  certain 
lifetime  and  give  off  photons  with  energy  near  the  band-gap  of  the  semiconductor. 
Nonradiative  processes  are  also  taking  place  that  allow  electron-hole  pairs  to  recombine 
without  producing  a  photon.  These  nonradiative  processes  reduce  the  number  of  carriers 
that  can  combine  radiatively,  and  generally  produce  additional  heat,  both  of  which  reduce 
the  overall  efficiency  of  opto-electric  devices.  The  nonradiative  recombination  methods 
that  will  be  discussed  are  Shockley-Read-Hall  (SRH)  and  Auger  recombination,  both  of 
which  are  capable  of  dominating  carrier  recombination.  The  recombination  rate  equation 
that  describes  how  the  carrier  density  will  change  over  time  can  be  written  as  (Agrawal 
andDutta,  1993): 

dn  2  3 

=  A  SRH  n+  B  rad-n  +  C  Auger'H 

^  ,  (2.5) 

where  n  is  the  carrier  density,  Asrh  is  the  SRH  coefficient.  Brad  is  the  radiative 
coefficient,  and  CAuger  is  the  Auger  coefficient. 

SRH  recombination  is  a  single  body  recombination  process  due  to  a  localized 
defect  that  creates  energy  levels  deep  within  the  band-gap  of  the  material  called  a  trap. 
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An  electron  can  give  up  its  energy  in  the  form  of  a  phonon  as  it  recombines  with  a  hole 
by  way  of  the  trap,  which  will  heat  the  lattice  and  will  not  produce  a  photon  (Coldren  and 
Corzine,  1995:144-148).  The  number  of  traps  in  a  semiconductor  will  vary  with  the 
impurities  and  lattice  vacancies  in  the  material.  SRH  recombination  will  be  the  dominant 
recombination  process  at  lower  carrier  densities. 

Unlike  SRH,  the  Auger  recombination  process  is  a  three-body  process,  and  will 
therefore  dominate  at  very  high  carrier  concentrations.  The  Auger  recombination  occurs 
when  the  energy  from  an  electron-hole  pair  recombination  is  transferred  to  another  carrier 
through  carrier-carrier  interactions  (Bhattacharya,  1997:147).  There  are  multiple 
combinations  in  which  Auger  recombination  can  occur,  however,  the  most  prominent  in 
this  material  is  called  Valence  Band  Auger  recombination  or  CHSH,  where  C  represents 
the  conduction  band,  H  the  heavyhole  band,  and  S  the  split-off  band  (Hausser,  1990). 

The  CHSH  Auger  process,  which  can  be  seen  in  the  schematic  in  Figure  2.3,  dominates 


Figure  2.3.  Valence  band  Auger  process  (CHSH).  An  electron-hole  pair  recombines 
and  the  resulting  energy  promotes  a  hole  from  the  split-off  band  to  the  heavy-hole  band. 
In  mid-IR  semiconductors.  Eg  ~  A  making  the  CHSH  process  dominate  at  high  carrier 
densities. 
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in  this  material  because  the  band-gap  energy  (Eg)  is  nearly  degenerate  with  the  energy 
difference  between  the  heavy -hole  and  split  off  bands  (A)  in  the  valence  band  structure. 

Significant  work  has  been  done  to  find  ways  to  suppress  Auger  recombination, 
which  dominates  at  high  carrier  concentrations.  Hausser  et  al.  compared  Auger 
recombination  rates  in  bulk  and  QW  InGaAs,  and  found  lower  Auger  coefficients  in  QW 
structures  (Hausser,  1990).  Strain  is  expected  to  affect  the  Auger  rate  because  it  can  alter 
the  energy  gaps.  Eg  and  A.  Auger  recombination  rates  were  found  to  vary  directly  with 
lattice  temperature  in  InAs/GalnSb  by  Jang  and  Elatte  (Jang  and  Elatte,  1998).  The 
energy  gaps.  Eg  and  A,  will  vary  with  temperature  differently,  causing  the  Auger  rate  to 
vary  with  lattice  temperature.  The  relationship  between  Auger  recombination  and  carrier 
temperature  will  be  discussed  in  Chapter  4.  It  is  evident  that  the  semiconductor  structure 
and  operating  temperatures  play  a  critical  role  in  the  Auger  process,  which  greatly  affect 
the  overall  efficiency  of  the  device. 

2,4,  Carrier  Scattering 

In  this  experiment,  a  pulsed  laser  is  used  to  create  electron-hole  pairs  so  that  the 
recombination  mechanisms  can  be  studied.  Electrons  are  excited  high  into  the 
conduction  band,  where  they  interact  with  each  other  and  the  lattice.  These  highly 
energetic  carrier  interactions  can  negatively  impact  radiative  recombination  rates  and 
overall  device  performance  (Jang  et  al,  1996).  The  dominant  interactions  that  will  be 
discussed  are  carrier-carrier  and  carrier-phonon  interactions. 

Carrier-carrier  scattering  allows  for  energy  exchange  among  the  highly  energetic 
carriers  created  with  an  excitation  pulse.  These  carriers  will  thermalize,  or  reach  a 
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common  temperature,  typieally  within  one  pieoseeond  of  the  exeitation  pulse  (Shah, 
1996:161).  Carrier-earrier  seattering  dominates  initially  but  eontributes  little  to  earner 
eooling  in  the  lattiee. 

The  main  souree  of  energy  transfer  from  the  earriers  to  the  lattiee  is  through 
earrier-phonon  seattering.  Carriers  will  lose  energy  to  phonons  in  the  lattiee,  eausing  the 
earner  temperature  to  eool  to  the  lattiee  temperature.  Previous  work  has  shown  that 
highly  energetie  earners  in  a  quantum  well  will  eool  through  interaetions  with  optieal 
phonons  at  earrier  temperatures  greater  than  40  K  (Yang,  1985).  A  MQW  strueture  made 
of  materials  similar  to  the  semieonduetors  in  this  experiment  was  found  to  eool  by  optieal 
phonon  emission  on  the  order  of  15  ps  (Jang  et  al,  1996).  The  eooling  of  energetie 
earriers  will  be  investigated  in  this  experiment  and  eompared  to  other  QW  materials.  The 
theory  of  carrier  dynamies  in  MQW  struetures  will  be  used  in  Chapter  4  to  analyze  data 
obtained  from  this  experiment,  and  ultimately  compare  performanee  eharaeteristies  of 
different  semieonduetor  samples. 
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3,  Experiment 


Luminescence  emitted  by  mid-IR  laser  devices  results  from  the  recombination  of 
electron-hole  pairs  excited  from  an  initial  laser  pulse.  Both  radiative  and  nonradiative 
recombination  of  electron-hole  pairs  takes  place  over  several  nanoseconds  after  the  initial 
excitation  pulse.  A  typical  method  for  collecting  time  resolved  spectra  uses  a  streak 
camera,  which  has  a  temporal  resolution  of  approximately  10  ps  (Shah,  1988).  In  order 
to  obtain  better  temporal  resolution,  Mahr  and  Hirsch  developed  the  upconversion 
technique  (Mahr  and  Hirsch,  1974).  Upconversion  is  applied  in  this  experiment  to  obtain 
sub-picosecond  resolution  of  the  luminescence  emitted  over  a  3-ns  time  period  (TRPL). 

A  mode  locked  Ti; Sapphire  laser  is  used  to  create  extremely  short  pulses.  The  beam  is 
split  and  one  leg  is  used  to  excite  a  semiconductor  sample.  The  luminescence  from  the 
excited  sample  is  focused  inside  a  nonlinear  crystal,  where  it  interacts  with  a  beam  that 
had  been  split  off  from  the  laser  source.  The  wave  mixing  that  occurs  within  the  crystal 
produces  a  third  beam,  with  a  predictable  wavelength,  that  is  collected.  The  portion  of 
the  beam  split  off  to  excite  the  sample  will  be  referred  to  as  the  PL  beam,  and  the 
luminescence  signal  will  be  called  the  PL  signal.  Together,  the  PL  beam  and  the  PL 
signal  travel  an  optical  path  that  is  equal  to  the  path  of  the  split  off  laser  beam,  which  will 
be  called  the  pump  beam.  The  third  beam,  called  the  upconverted  beam,  will  only  be 
created  when  both  the  PL  signal  and  pump  beam  are  present  within  the  crystal.  A  pulsed 
pump  beam  will  therefore  act  as  a  light  gate  and  allow  for  the  temporal  resolution  of  the 
PL  signal  (Shah,  1988).  By  varying  the  optical  path  of  one  of  the  incoming  beams,  the 
luminescence  from  the  excited  semiconductor  can  be  resolved  over  time.  This  process  is 
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depicted  in  Figure  3.1.  This  chapter  will  discuss  the  theory  of  upconversion,  the  methods 
applied  in  the  experiment  to  gather  the  data,  and  the  efficiency  of  the  upconversion 
process.  It  will  conclude  with  a  description  of  the  sample  semiconductors  investigated. 


Figure  3.1.  Wave  mixing  resulting  in  the  time  resolution  of  the  luminescence. 

(a)  Depicts  the  luminescence  curve  over  time,  (b)  Represents  the  pump  beam  at 
some  delay  time  x. 

3,1.  Upconversion  theory 

Upconversion,  also  called  sum  frequency  generation,  is  the  process  of  mixing  two 
optical  waves  to  generate  a  third  wave.  The  interaction  required  for  the  upconversion  of 
two  optical  waves  occurs  in  a  non-linear  optical  medium.  A  non-linear  medium  is  one  in 
which  the  polarization  density  can  be  expressed  as  a  Taylor  series  expansion: 

Pi  =  s  Xij  Ej  +  2s  x,ijk  E,j  Ek  +  . . .,  (3.1) 

where  P  is  the  polarization  density,  E  is  the  induced  electric  field,  s  is  the  permeability  of 
free  space,  andx  is  the  susceptibility  (Yariv,  1984:504).  In  this  experiment,  waves  are 
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mixed  in  a  KTi0As04(KTA)  crystal  exhibiting  second-order  nonlinearity,  and  higher 
order  terms  in  the  series  expansion  are  considered  negligible.  If  the  optical  waves  are 
modeled  as  plane  waves,  the  second  order  term  in  Eq  (3.1)  results  in  two  conditions 
required  for  upconversion.  The  first  is  the  frequency  matching  condition  that  dictates: 

(03  =  C0i  +  C02  (3^ 

In  this  experiment,  the  pump  beam  frequency  and  the  peak  luminescence 
frequency  from  the  excited  sample  sum  to  the  upconverted  frequency.  Manipulating  Eq 
(3.2)  creates  an  expression  based  on  the  conservation  of  energy  used  to  determine  the 
wavelength  of  the  upconverted  photons: 

_[ _ \_  1 

^  up  ^  PL  ^  Pump  /3  - 


where  kup  is  the  wavelength  of  the  upconverted  wave  and  kpL,  and  kpump  refer  to  the 
wavelengths  of  the  PE  and  pump  beams,  respectively. 

The  second  condition  required  for  upconversion  is  the  phase  matching  condition: 

k  up  =  k  PL  -I-  k  Pump  /3  a 


where  k  is  the  wave  vector  ( |k|  =  27in  /  X)  illustrated  in  Eigure  3.2.  When  the  law  of 
cosines  is  applied  to  Eq  (3.4)  and  k  is  replaced,  the  expression  becomes: 


Hup 
^  up 


2  2 

npL  npump  Hup-npL 
- 1 - 2 - cos 


I 


PL 


I 


Pump 


^  up-^  PL 


(ti  -  y) 


(3.5) 


where  y  is  the  angle  between  the  PL  and  pump  wave  vectors  and  n^p,  npL,  and  np„mp  are 
the  refraction  indices  seen  by  the  upconverted  beam,  the  PL  signal,  and  the  pump  beam. 
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Figure  3.2.  Phase  matching  geometry  required  for  upconversion.  The  pump 
beam  and  PL  beam  enter  the  crystal  at  20° 


Fenimore  (1995)  gives  the  Sellmeier  equations,  which  express  index  of  refraction  as  a 
function  of  wavelength  for  the  KTA  crystal  in  terms  of  the  three  principal  axes  nx,  ny,  and 
iiz.  The  crystal  used  here  is  cut  at  an  angle  O  =  43°,  where  O  is  the  angle  between  the 
crystal  axis  and  the  normal  to  the  surface  in  the  X-Z  plane.  The  refraction  indices  in  the 
ordinary  and  extraordinary  polarizations  can  now  be  written  in  terms  of  the  crystal  axis 
indices: 

noU)  =  ny(z) 


n 


nx(x)-ny(z) 

•cos(0)^+  nx(^)^-sin(0)^ 


(3.6) 


where  0  is  the  angle  between  the  extraordinary  beam  and  the  z  axis.  When  the  laser  beam 
is  split,  the  PL  beam  is  orientated  to  be  ordinary  and  the  pump  beam  is  extraordinary. 

The  properties  of  the  KTA  crystal  dictate  that  the  upconverted  wave  will  then  be 
ordinary.  With  all  other  parameters  known,  Eq  (3.6)  can  be  solved  for  0  for  any  given 
wavelength.  Recall  that  0  is  the  angle  of  the  pump  beam  to  the  z-axis  and  can  be  changed 
by  simply  rotating  the  crystal.  A  MathCad  program  (See  Appendix  A.l)  incorporates 
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these  equations  so  that  it  will  provide  the  erystal  tilt  required  to  meet  the  phase  matehing 
eondition  for  a  given  experimental  set  up.  The  next  seetion  will  deseribe  how  the 
experiment  was  eonstrueted  to  mix  these  waves  inside  a  erystal  and  eolleet  the  resulting 
upeonverted  beam. 

3,2,  Experimental  Set  Up 

Capt  William  Cooley  first  developed  the  upeonversion  experiment  at  AFIT  in 
1997.  A  sehematie  of  the  experimental  set  up  ean  be  seen  in  Figure  3.3.  This  seetion 
details  the  laser  souree  used  for  exeitation,  the  optieal  path,  and  the  eolleetion  of  the 
photolumineseenee  and  upeonverted  signal. 

3,2,1,  Laser  Source 

A  Coherent  Mira  900  mode-loeked  TiiSapphire  laser  operating  at  1.5  watts 
average  power  produeed  the  pulses  required  for  upeonversion.  The  laser’s  repetition  rate 
is  76  MHz  with  a  eenter  wavelength  of  809  nm.  The  mode-loeked  pulses  exeite  the 
sample  every  13.2  ns,  allowing  time  to  study  the  lumineseenee  deeay.  The  pulse  width, 
whieh  establishes  the  temporal  resolution  of  the  lumineseenee  deeay,  is  122  fs.  The 
Ti:Sapphire  is  optieally  pumped  by  a  12.0  watt  Speetra  Physies  argon-ion  (Ar^)  laser. 

The  speetral  and  temporal  aspeets  of  the  Ti:Sapphire  beam  ean  be  seen  in  Figure  3.4. 

The  pulse  width  was  measured  using  an  INRAD  Model  5-14BX  Autoeorrelator 
and  verified  using  the  upeonversion  experiment.  The  autoeorrelator  eould  produee  highly 
aeeurate  results  mueh  faster  than  the  upeonversion  experiment.  The  autoeorrelator  splits 
the  laser  into  two  paths  and  foeuses  the  beams  together  inside  a  non-linear  erystal. 
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Figure  3.3.  Experimental  set  up  in  TREE  experiment.  The  Ti; Sapphire  laser 
ereates  a  122  fs  pulse  that  is  split  into  two  paths.  The  PE  beam  travels  through  a 
delay  stage  before  exeiting  the  sample.  Eumineseenee  is  eolleeted  and  mixed  with 
the  pump  beam  in  the  KTA  erystal.  The  upeonverted  beam  is  then  direeted  to  the 
speetrometer  and  photon  eounter. 
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A  LilOa  crystal  is  used  to  create  an  upeonverted  beam.  One  path  passes  through  a  delay 
stage  that  varies  the  optical  path  length,  while  a  detector  collects  the  upeonverted  beam. 
The  result,  whieh  ean  be  seen  in  Figure  3.4,  is  a  temporal  pulse  with  a  measurable  width. 
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Figure  3.4.  (a)  Temporal  pulse  width  of  Ti:Sapphire  laser.  FWHM  is  measured  to 
be  Atm  =188  fs.  A  correction  factor  of  0.648  is  used  to  obtain  an  actual  pulse  width 
of  Ata  =  122  fs  (b)  Spectral  width  of  TiiSapphire  pulse  with  a  FWHM  of  10  nm. 
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The  manufacture  suggests  using  a  sech  function  to  fit  the  line  shape,  which  gives  rise  to 
a  conversion  factor  of  0.648  to  convert  from  the  measured  pulse  width  to  the  actual  pulse 
width. 

Knowledge  of  the  focused  laser  spot  size,  both  within  the  crystal  and  on  the 
sample,  is  critical  for  performing  the  upconversion  experiment  and  is  used  extensively  in 
the  data  analysis  seen  in  Chapter  4.  A  Coherent-Mode  Master  was  used  to  measure  the 
beam  diameter  as  well  as  the  M  value  of  the  beam.  The  M  value  is  a  relative  term 
comparing  the  beam  profde  to  that  of  a  Gaussian  profde.  The  Ti: Sapphire  laser  produced 
a  beam  with  an  M  value  of  1 .3,  which  means  the  beam  waist  will  be  1 .3  times  the  size  of 
a  diffraction  limited  Gaussian  beam  waist.  The  beam  diameter  was  measured  in  two 
places  in  the  beam  path  to  be  6.85  +/-  .032  mm.  The  lens  used  to  focus  the  pump  beam 
and  PL  beam  on  the  crystal  have  a  focal  length  of  200  mm  and  300  mm,  respectively. 
Assuming  the  beam  is  Gaussian,  the  spot  sizes  on  the  crystal  can  be  calculated  by: 

4  ?i-f 

2(0  = - 

^  ^  ,  (3.7) 

where  2co  is  the  focused  beam  diameter,  f  is  the  focal  length,  and  D  is  the  beam  diameter 
just  before  reaching  the  lens.  Using  this  method,  the  pump  beam  spot  size  was 
determined  to  be  39.1  +/-  0.15  pm  and  the  PL  spot  size  was  58.5  +/-  0.20  pm.  The  M^ 
value  of  the  beam  was  incorporated  in  the  spot  size  calculations  to  account  for  the  beam’s 
shape.  An  analysis  of  the  accuracy  required  of  the  spot  size  measurement  is  included  in 
Chapter  4.  When  operating  properly,  the  TiiSapphire  laser  produced  a  very  stable  beam 
with  consistent  pulse,  which  proved  crucial  in  day-to-day  data  collection. 
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3.2.2,  Optical  Path 

Once  the  beam  leaves  the  Ti:Sapphire  laser,  wavelengths  below  700  nm  are 
filtered  out,  and  the  beam  is  expanded  and  collimated  before  it  reaches  a  beam  splitter,  as 
seen  in  Figure  3.3.  A  7,/2  wave  plate  is  used  to  vary  the  polarization  of  the  beam.  The 
beam  splitter  separates  the  beam  into  horizontal  and  vertical  components.  By  adjusting 
the  yi  wave  plate,  the  power  ratio  in  eaeh  leg  ean  be  varied.  The  PL  beam  is  the  portion 
that  will  excite  the  sample  creating  a  PL  signal.  The  pump  beam  will  propagate  a  fixed 
distance  eomparable  to  the  combined  distance  the  PL  beam  and  PL  signal  will  travel. 

The  PL  signal  and  pump  beam  will  be  focused  onto  the  nonlinear  crystal  causing  them  to 
mix  and  produee  an  upeonverted  beam. 

3.2.2,1.  PL  Path 

The  PL  path  is  eomprised  of  the  PL  beam  and  the  PL  signal.  The  PL  beam  begins 
when  the  vertically  polarized  beam  exits  the  beam  splitter,  passes  through  a  neutral 
density  fdter,  and  eneounters  a  delay  stage  seen  in  Figure  3.3.  The  delay  stage  eonsists  of 
a  gold-coated  retrorefieetor  on  a  sliding  mount  that  allows  the  optieal  path  length  of  this 
leg  to  be  varied.  The  ability  to  change  the  optieal  path  of  one  beam  with  respeet  to  the 
other  allows  the  luminescence  decay  to  be  mapped,  as  described  by  Figure  3.1.  The 
beam  makes  two  passes  through  the  30-em-long  delay  stage,  whieh  allows  the  path  length 
of  the  PL  beam  to  vary  by  4  ns  (4  x  30  em  /  c  =  4  ns).  The  PL  beam  is  then  focused  onto 
the  sample  being  investigated,  which  is  mounted  in  a  vaeuum  ehamber.  A  gold-eoated, 
90°  off-axis  parabolic  mirror  with  a  focal  length  of  150  mm  directs  the  luminescence,  or 
PL  signal,  from  the  sample  to  an  identical  mirror,  which  focuses  the  luminescence  on  the 
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nonlinear  crystal  mounted  on  the  rotation  stage.  The  two  mirrors  form  a  one-to-one 
imaging  system,  with  the  sample  in  the  object  plane  and  the  crystal  in  the  image  plane. 

3.2.2.2,  Pump  Beam  Path 

The  pump  beam  is  horizontally  polarized  when  it  exits  the  beam  splitter  and 
travels  over  an  optical  path  similar  in  length  to  the  PL  path.  A  200-mm  lens  is  used  to 
focus  the  beam  down  to  a  spot  within  the  nonlinear  crystal,  as  seen  in  Figure  3.3.  A 
prism  is  used  to  create  a  20°  angle  between  the  pump  beam  and  the  PL  beam.  The 
relative  angle  between  the  beams  is  fixed  and  used  in  Eq  (3.5)  when  determining  the 
crystal  tilt  angle  required  for  upconversion.  The  pump  beam  is  mixed  with  the  PL  signal 
to  create  the  upconverted  beam. 

3.2.2.3,  Beam  Overlap 

The  upconversion  process  relies  on  two  optical  waves  interacting  with  one 
another  to  generate  a  third  wave  (Saleh  and  Teich,  1991:748).  This  requires  that  the 
waves  overlap  in  space  and  time  in  order  to  create  an  upconverted  beam.  Spatial  overlap 
was  achieved  by  replacing  the  non-linear  crystal  with  a  pinhole  in  a  similar  mount  on  the 
rotation  stage  seen  in  Figure  3.3.  Both  the  pump  beam  and  the  PL  signal  were  focused 
through  a  100-pm  pinhole  initially  to  determine  coarse  alignment.  Although  the  beam 
spot  sizes  were  calculated  to  be  39.1  and  58.5  pm  for  the  pump  beam  and  PL  signal,  the 
100-pm  pinhole  did  not  ensure  sufficient  overlap  for  upconversion.  A  25-pm  pinhole 
was  then  used  to  make  fine  adjustments  and  ensure  maximum  spatial  overlap  of  the  two 
beams. 
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Temporal  overlap  was  aehieved  by  mixing  the  pump  beam  with  itself  to  ereate  an 
upeonverted  beam.  The  sample  was  moved  out  of  the  PL  beam  path,  and  the  laser’s 
specular  reflection  off  of  the  cold  linger  in  the  sample  chamber  was  collected  rather  than 
the  luminescence.  This  process  is  identical  to  the  method  used  by  the  autocorrelator 
discussed  earlier,  and  a  LilOs  crystal  is  used  to  mix  the  two  beams.  The  resulting 
upeonverted  signal  is  in  the  visible  spectrum,  so  that  it  is  easily  seen  when  the  optical 
path  lengths  of  the  two  legs  are  equal.  The  delay  stage  was  moved  incrementally  until 
the  upeonverted  beam  was  observed.  This  procedure  confirmed  that  both  beams  were 
present  in  the  same  space  at  the  same  time  in  the  crystal. 

3.2,4.  Signal  Collection 

This  experiment  requires  that  standard  luminescence  data  be  collected  on  each 
sample  to  identify  the  spectral  characteristics  of  the  sample’s  emission  prior  to 
performing  upconversion.  The  design  of  the  experiment  allows  for  luminescence  data  to 
be  taken  with  minimal  impact  to  the  upconversion  configuration. 

3,2,4.1.  Photoluminescence  experiment 

No  wave  mixing  is  required  when  taking  time-averaged  photoluminescence  data. 
The  nonlinear  crystal  is  removed  from  the  rotation  stage  in  Figure  3.3  and  the  pump  beam 
is  blocked.  The  PL  beam  excites  the  sample  and  parabolic  mirrors  direct  the 
luminescence.  The  PL  signal  is  then  collected  and  re-collimated  by  a  CaF2  lens  that  will 
pass  light  in  the  IR  spectrum.  A  second  CaF2  lens  focuses  the  signal  on  the  entrance  slit 
of  a  V2  meter  SPEX  spectrometer  with  a  4-pm-blaze  grating.  The  grating  has  1200 
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grooves/mm  with  a  calculated  resolution  of  5  nm.  A  Standard  Research  Lock-in 
Amplifier  amplifies  the  signal  from  a  single  element  InSb  detector. 

3.2.4,2,  TRPL  experiment 

The  nonlinear  crystal  used  to  mix  the  PL  and  pump  beams  in  this  experiment  is  a 
1-mm-thick  KTA  crystal  cut  so  the  z-axis  is  43°  from  the  incident  normal.  The  pump  and 
PL  beams  must  overlap  in  time  and  space  within  the  erystal  for  them  to  mix  and  create 
the  upconverted  beam.  With  all  other  parameters  fixed,  the  wavelength  of  the 
luminescenee  signal,  typically  3.1-  4.5  pm,  dictates  the  erystal  tilt  angle  and  the  exit 
angle  of  the  upconverted  beam.  A  collection  lens  is  used  to  focus  the  upconverted  signal 
onto  the  entrance  slit  of  a  14  meter  SPEX  speetrometer.  The  grating  used  has  1200 
grooves/mm  with  a  calculated  resolution  of  3.3  nm.  The  signal  is  collected  by  a  GaAs 
photomultiplier  tube  cooled  to  0°  C  and  biased  at  2000  V.  The  signal  is  amplified  and 
sent  to  a  Standard  Research  SR400  photon  counter,  which  passes  data  to  a  computer 
interface  used  to  control  the  experiment.  The  computer  allows  for  either  the  speetrometer 
or  the  delay  stage  to  be  incremented  while  taking  photon  counts. 

In  order  to  reeord  the  luminescenee  temporal  evolution,  the  speetrometer  was 
tuned  to  the  peak  upconverted  wavelength,  and  the  crystal  angle  was  set  to  maximize  the 
upconverted  signal.  The  delay  stage  was  then  moved  inerementally  as  intensity  data  was 
recorded.  The  photon  counter  integrated  the  counts  it  received  over  30  seconds  at  each 
increment  in  the  delay.  When  investigating  the  decay  of  the  luminescence,  the 
experiment  was  set  to  record  data  every  50  ps  over  a  3  ns  time  period.  The  number  of 
data  points  allowed  for  a  high  confidence  eurve  fit  to  analyze  the  luminescence  decay. 
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More  data  points  were  required  to  investigate  the  lumineseenee  rise,  whieh  oeeurs  over 
the  first  25-30  ps.  Data  was  reeorded  every  2  ps  over  a  50  ps  time  period  to  examine  the 
rise  time. 

Speetral  data  was  taken  by  setting  the  delay  stage  at  one  position  and  seanning 
over  the  speetra  of  the  upeonverted  signal.  The  speetrometer  sean  speed  was  dietated  by 
the  resolution  of  the  grating  ealeulated  to  be  3.3  nm.  In  order  to  mateh  the  grating 
resolution,  a  speed  of  10  nm/min  was  ehosen  so  that  eaeh  data  point  would  represent  the 
intensity  over  5  nm.  This  proeedure  was  repeated  for  several  fixed  times  over  the  3 -ns 
period  after  the  exeitation  pulse. 

Another  method  to  obtain  speetral  data  was  also  attempted  that  has  been  used  in 
previous  experiments  (Franz,  1997:33).  The  upeonversion  experiment  was  run  to 
measure  the  time  deeay  of  the  lumineseenee  at  a  fixed  wavelength.  The  wavelength  was 
then  ineremented  and  another  time  deeay  plot  was  taken.  All  of  the  points  eorresponding 
to  a  speeifie  delay  time  were  then  put  together  to  ereate  a  speetral  plot  at  a  snapshot  in 
time.  This  method  produees  signilioantly  more  data  points  in  time,  but  many 
measurements  would  have  to  be  taken  to  obtain  a  suffioient  number  of  speetral  points. 
The  analysis  of  this  data  focuses  on  the  spectral  data,  which  is  improved  with  more  data 
points.  In  addition,  the  phase  matching  condition  to  upconvert  specific  wavelengths 
required  a  crystal  tilt  angle  unavailable  in  this  experimental  set  up.  If  the  crystal  is  tilted 
between  18.8°  and  20.5°,  the  specular  reflection  retraces  the  path  of  the  pump  beam  and 
interferes  with  the  mode-locking  of  the  TiiSapphire. 
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3,3.  Upconversion  Efficiency 

When  the  upeonversion  experiment  is  used  to  measure  the  time  deeay 
luminescence,  the  crystal  is  turned  to  the  phasematching  angle  in  order  to  satisfy  Eq 
(3.4),  which  will  maximize  upconversion.  As  the  crystal  is  titled  away  from  the 
phasematching  angle  so  that  Eq  (3.4)  is  no  longer  true,  the  upconversion  efficiency 
decreases  rapidly  until  upconversion  is  no  longer  possible.  Spectral  data  was  taken  in 
this  experiment  by  fixing  the  crystal  to  the  optimal  phasematching  angle  to  upconvert  the 
peak  PE  signal.  In  order  to  ensure  that  the  spectral  tails  of  the  luminescence  can  be 
upconverted,  the  acceptance  angle  of  the  crystal  must  be  calculated.  The  acceptance 
angle  is  the  solid  angle,  which  is  centered  about  the  phasematching  angle,  over  which 
upconversion  can  still  take  place.  The  acceptance  angle  is  given  by: 


= 


2.78-n-A-  pi 


(3.8) 


where  E  is  the  crystal  thickness  (Shah,  1988).  The  index  of  refraction  for  the  PE  signal 
(n  =  1.74)  and  the  upconverted  signal  (nup=1.85)  are  calculated  using  the  program  in 
Appendix  A.l.  The  sample  with  the  broadest  spectra,  spanning  from  3.1  to  3.7  pm,  was 
chosen  to  analyze  the  acceptance  angle.  The  calculated  acceptance  angle  for  Sample  91 
is  A(p  =  0.021  steradians.  This  solid  angle  is  converted  to  a  planar  angle  using  the 
following  equation: 

4)  =  2-7i(i  -  cos  (e  max))  /o  ( 
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The  resulting  planar  angle  that  will  be  aecepted  by  the  crystal  for  upconversion  is 
+/-  4.6°.  The  phasematching  angles  for  the  Sample  91  PL  spectra  are  19.1°  +/-  3.5°, 
which  is  within  the  acceptance  angle  for  all  wavelengths  in  the  spectra.  While  the  entire 
spectra  can  be  upconverted,  the  efficiency  of  the  upconversion  process  will  decrease  as 
the  crystal  tilt  angle  moves  away  from  the  phasematching  angle. 

Cooley  calculated  the  upconversion  efficiency  for  this  experiment  to  be 
approximately  0.1%  for  perfect  phasematching  of  the  incoming  beams  (Cooley, 

1997:81).  This  peak  efficiency  is  based  on  the  crystal  parameters,  the  wavelengths  being 
mixed,  and  the  ratio  of  beam  power  to  area  on  the  crystal,  and  should  remain  constant  for 
a  specific  experimental  set  up.  When  the  phasematching  condition  is  no  longer  met,  Eq 
(3.4)  becomes: 


Ak  =  k  PL  +  k  Pump  -  k 


up 


(3.10) 


The  upconversion  efficiency  will  fall  off  as  (Shah,  1988): 

ri(Ak)  _  sin(L  g-Ak)^ 

^  peak  (ke-Ak)^ 

9 


where  Eg  is  the  effective  length  of  the  crystal,  which  may  be  shorter  than  the  actual  length 
of  the  crystal. 

Dmitriev’s  method  for  determining  the  effective  length  for  an  upconversion 
experiment  with  short  laser  pulses  is  to  compare  the  crystal  length  to  the  quasi-static 
length  given  by: 


E 


qsp 


(v-Aco)  ^ 


(3.12) 
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where  Aco  is  the  speetral  pulse  width  and  v  is  the  inverse  group  velocity  mismatch  given 
by: 


-  1  - 1 
V  =  u 1  -U2 


(3.13) 


where  ui  and  U2  are  the  group  velocities  for  the  two  waves  being  mixed  in  the  crystal 
(Dmitriev  et  al.,  1997:60).  The  calculated  quasi-static  length  was  1 12  pm,  which  is 
significantly  shorter  than  the  actual  crystal  length  of  1  mm.  Cooley  calculated  the 
effective  length  for  this  experiment  to  be  Le  =  106  pm,  taking  into  account  all  spatial  and 
temporal  considerations  (Cooley,  1997:82-84).  The  effective  length  used  will  be  the 
quasi-static  length  (Le  =  112  pm)  because  the  simplified  approach  proposed  by  Dmitriev 
produces  results  agreeable  with  Cooley. 

The  effective  length  can  now  be  used  in  Eq  (3.1 1)  to  determine  the  upconversion 
efficiency  for  a  particular  experiment.  The  sample  with  the  widest  spectral  width  was 
used  to  calculate  Ak  values  for  various  wavelengths  within  the  PL  spectra.  The 
relationship  between  the  upconversion  efficiency  fall  off  and  wavelength  for  a  particular 
semiconductor  can  be  seen  in  Figure  3.5.  The  upconversion  efficiency  is  greater  than 
97%  at  the  FWHM  of  the  spectrum  and  greater  than  90%  across  the  entire  spectrum.  The 
variation  in  efficiency  will  introduce  error  in  spectral  luminescence  data,  but  is  small 
enough  to  allow  for  qualitative  comparisons  of  carrier  temperature  derived  from  this  data. 


3,4,  Sample  Description 

The  four  samples  investigated  in  this  experiment  were  laser  structures  provide  by 
Lincoln  Laboratory  at  the  Massachusetts  Institute  of  Technology  (MIT).  Specific 
information  on  growth  procedures  has  been  published  (Turner,  1995).  Cooley  (1997)  and 
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McKay  (2000)  have  each  investigated  some  of  the  samples  tested  in  this  experiment.  The 
rationale  behind  examining  these  semiconductors  is  to  provide  a  uniform  comparison  of 
different  type  structures  using  a  consistent  method.  In  addition,  data  involving 
coefficient  reliability,  luminescence  rise  time,  and  spectral  decay  was  collected  in  order 
to  better  understand  the  carrier  dynamics  in  each  sample. 


Figure  3.5.  Upconversion  efficiency  compared  to  peak  efficiency  over  the  luminescence 
spectra  of  Sample  91.  The  Sample  91  PL  curve  peaks  at  3.5  pms  with  a  FWHM 
indicated  by  the  dashed  lines  in  the  figure. 

Sample  “B”  is  a  multiple  quantum  well  type  I  structure  grown  by  molecular  beam 
epitaxy  (MBE).  The  confinement  of  electrons  and  holes  occurs  in  the  same  material  in  a 
type  I  structure.  The  well  structure  is  made  up  of  InAso.gssSbo.oes/Ino.gsAlo.isAso.gSbo.i  and 
can  be  seen  in  Figure  3.6. 
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Figure  3.6.  Sample  B  quantum  well  strueture.  Energy  diagram  with  offset  bands  and 
exeitation  energy  drawn  to  seale.  Well  and  Barrier  widths  are  also  ineluded. 

Samples  88,  90,  and  91  are  InAs/Ga.gsIn.isSb  type  II  quantum  well  struetures  also 
provided  by  MIT.  A  Type  II  strueture  allows  for  eleetron  eonfinement  in  one  material, 
while  holes  are  eonfined  in  the  neighboring  material.  The  different  samples  are  identieal 
to  eaeh  other  exeept  for  varying  well  widths,  whieh  ean  be  seen  in  Figure  3.7.  Time 
integrated  photolumineseenee  was  taken  on  eaeh  of  the  samples  and  ean  be  seen  in  Figure 
3.8.  The  data  was  taken  with  the  sample  held  at  77  K  and  an  exeitation  power  of  600 
mW.  Sample  B  and  sample  91  were  the  brightest  samples.  As  the  width  of  the  wells  in 
the  type  II  strueture  beeomes  wider,  the  peak  wavelength  shifts  to  lower  energies.  This 
result  is  predieted  by  Eq  (2.1),  whieh  states  that  energy  levels  in  a  quantum  well  are 
inversely  proportional  to  well  width. 


32 


Ga.gsIn.isSb 


Sample 

Well  Width 

88 

26  A 

90 

23  A 

91 

20  A 

24  A 


20  wells 


Figure  3.7.  Sample  88,  90,  and  91  quantum  well  structure.  Energy  diagram  with  offset 
bands  and  excitation  energy  drawn  to  scale.  Well  and  Barrier  widths  are  also  included. 
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Figure  3.8.  Photoluminescence  curves  for  each  sample  tested.  Sample  B  is  a  type  I 
structure  Sample  88,  90,  and  91  are  type  II.  Intensities  are  shown  as  true  scale. 
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4,  Data  Analysis  and  Results 


The  data  collected  in  the  upconversion  experiment  described  in  the  previous 
chapter  provides  both  time  resolved  luminescence  and  spectral  luminescence  intensities 
radiated  by  the  tested  samples.  Time  resolved  plots  are  used  to  determine  the 
recombination  coefficients  and  allow  for  a  qualitative  comparison  of  carrier  relaxation 
into  the  quantum  well  structure.  Spectral  data  is  used  to  extract  carrier  temperature 
information.  Recombination  rates,  carrier  relaxation,  and  carrier  temperature  are  all 
critical  factors  that  influence  laser  performance.  This  chapter  will  discuss  the  time 
resolved  and  spectral  results,  and  the  analysis  performed  on  the  experimental  data  will  be 
used  to  draw  conclusions  on  the  device  performance. 

4,1.  Determining  Recombination  Coefficients 

The  time  resolved  luminescence  decay  curves  recorded  over  3-ns  are  used  to 
determine  the  recombination  coefficients  for  the  samples.  As  mentioned  earlier,  the 
Auger  coefficient  is  of  greatest  importance  because  it  is  detrimental  to  device 
performance  and  can  be  reduced  by  altering  the  design  of  the  structure  (Lindle,  1995). 
There  have  been  some  discrepancies  in  calculating  recombination  rates  from 
luminescence  data  in  previous  reporting,  with  Auger  rates  ranging  from  10'  cm  /s 
(Lindle,  1995)  to  10'  cm  /s  (Cooley,  1997)  for  mid-IR  semiconductor  type  I  structures 
and  from  10'^^  cm%  (Meyer,  1998)  to  10'^^  cm%  (McKay,  2000)  in  type  II  structures. 
The  factors  that  lead  to  the  greatest  uncertainty  in  calculating  recombination  coefficients 
when  conducting  time  resolved  luminescence  experiments  arise  when  determining  carrier 
densities  and  curve  fitting  the  data.  Approximations  to  unknown  material  parameters 
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have  a  significant  impact  on  calculating  carrier  densities,  and  result  in  recombination 
coefficients  from  different  experiments  that  cannot  be  compared.  The  theory  remains 
sound  and  allows  for  different  materials  to  be  compared  provided  the  calculations  are 
subject  to  identical  assumptions  and  error  is  minimized.  This  section  deals  with  the 
analysis  performed  to  determine  recombination  coefficients  and  the  error  associated  with 
assumptions  that  have  been  made. 


4,1.1,  Analysis  of  Decay  Curve 

The  equation  that  describes  the  change  in  carrier  density  due  to  spontaneous 
emission  in  a  semiconductor  is: 

dn  .  „  2  „  3 

—  ^  SRH'm"  1^  rad ^  Auger 

^  ,  (4.1) 

where  n  is  the  carrier  density,  and  Asrh,  Brad,  and  CAuger  are  coefficients  associated  with 
the  recombination  processes  discussed  in  Chapter  2.  Eq  (4.1)  can  be  solved  for  carrier 
density  as  a  function  of  time  if  the  initial  condition,  or  the  carrier  density  at  time  t  =  0,  is 
known. 


The  data  collected  in  the  upconversion  experiment  yields  luminous  intensity  at 
various  times,  as  depicted  in  Figure  4.1. 

A  relationship  between  carrier  density  and  luminescence  allows  the  collected  data 
to  be  transformed  to  carrier  densities  at  the  same  time  intervals.  Eq  (4.1)  can  then  be  fit 
to  the  data,  and  the  coefficients  extracted.  These  coefficients  are  parameters  of  the 
material  and  allow  for  comparison  between  different  quantum  well  types  and  sizes. 

The  carrier  densities  must  be  known  before  a  relationship  between  carrier  density  and 
luminescence  can  be  determined.  The  method  employed,  which  is  taken  from  Sermage 
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Delay  (psec) 

X  300  mW 
-H  225  mW 
n  150  mW 
O  75  mW 


Figure  4.1.  Sample  TRPL  data  collected.  The  intensity  of  the  upconverted  signal 
from  an  optically  excited  quantum  well  structure  is  plotted  over  3  nanoseconds. 
Multiple  curves  indicate  varying  excitation  power. 


(Sermage,  1997)  and  adapted  by  Cooley,  allows  the  maximum  carrier  density  created  in 
the  structure  to  be  written  as  a  function  of  laser  input  power,  beam  sizes,  and  material 
parameters: 


E 

ht) 


^Nqw'Tt'Wpl  -tej 


®  'N  qw'(t  B+t  bar) 


1 


2 

Wp 

1  +  — ^ 
2 

Wpl 


(4.2) 


The  initial  term  represents  the  number  of  photons  per  pulse.  E  is  the  energy  per 
pulse  and  is  determined  by  taking  the  average  power  divided  by  the  repetition  rate  of  the 
laser.  The  second  term  represents  the  number  of  carriers  per  unit  volume.  Nqw  is  the 
number  of  quantum  wells  in  the  sample,  te  is  the  thickness  of  the  quantum  well,  and  Wpi 
and  Wp  are  the  waist  sizes  of  the  PL  beam  and  pump  beam  defined  by  the  1/e  intensity 
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points.  The  third  term  indicates  the  power  loss  associated  with  reflection  off  the  sample. 
R  is  the  reflectance  and  can  be  measured  directly  or  calculated  by: 


Rs 


rns-iY 

^ns+  1  j 

9 


(4.3) 


where  Ug  is  the  index  of  refraction  of  the  semiconductor.  The  fourth  term  represents  the 
amount  of  energy  absorbed  by  each  quantum  well.  The  term  a  represents  the  absorption 
coefficient  and  tbaris  the  barrier  layer  thickness.  The  final  term  is  a  beam  overlap 
efficiency  required  because  the  PL  and  pump  beam  are  not  focused  to  the  same  size, 
resulting  in  a  loss  of  photons  available  for  upconversion.  This  equation  was  modified  for 
type  II  structures  (McKay,  2000)  to  account  for  the  offset  wells: 


n(to)  = 


h-n 


t  InAs  GalnSb 


qw'Tt'Wpl  't  inAs't  GalnSb  ) 


(l-Rs)-i-e 


■  a  •  N  qYv  •  (t  InAs  GalnSb) 


W  1 


1  + 


w 


pi 


(4.4) 


where  tinAs  and  toainSb  are  the  thickness’  of  corresponding  layers  in  the  material. 

The  initial  carrier  density  equation  is  a  large  contributor  to  the  discrepancies  in 
reporting  recombination  coefficients  on  specific  materials.  The  initial  carrier  density 
equations  described  in  Eq  (4.2)  and  Eq  (4.4)  make  some  key  assumptions.  As  Cooley 
pointed  out,  the  carrier  loss  that  occurs  during  the  pulse  is  considered  negligible,  and 
hence  ignored  (Cooley,  1997:1 17).  The  closest  index  of  refraction  for  the  substrate  GaSb 
is  determined  experimentally  and  is  an  approximation  to  the  index  of  the  QW  structure 
(Madelung,  1991).  The  beam  overlap  was  calculated  assuming  a  Gaussian  beam  profile. 
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Finally,  the  absorption  coefficient  is  determined  by  ellipsometry  for  Sample  B,  and 
unavailable  for  Sample  88,  90,  and  91.  The  uncertainties  in  these  values  lead  to 
uncertainties  in  the  conversion  of  luminescence  to  carrier  densities,  and  therefore,  the 
recombination  coefficients. 

The  error  in  estimating  the  carrier  density  is  investigated  by  determining  the 
impact  each  parameter  in  Eq  (4.2)  has  on  the  computed  initial  carrier  density.  There  are 
uncertainties  in  the  absorption  coefficient,  index  of  refraction,  and  beam  waist  size  that 
have  been  decreased  through  experimentation  and  a  review  of  the  literature.  A  range  of 
absorption  coefficients  was  used  to  determine  carrier  densities,  and  the  percent  change  in 
carrier  density  over  that  range  was  calculated.  This  process  was  repeated  for  each 
unknown  parameter,  and  the  results  can  be  seen  in  Table  4.1.  This  analysis  indicated  that 
the  uncertainty  in  the  absorption  coefficient  has  the  strongest  influence  on  the  carrier 
density. 

While  the  absorption  coefficient  was  determined  experimentally  for  Sample  B,  no 
information  is  available  for  type  II  structures.  An  absorption  coefficient  of  2  x  10^  cm'^ 
was  determined  by  ellipsometry  on  Sample  B  (Cooley,  1997:120).  This  value  was  used 
in  the  analysis  for  all  four  samples.  It  is  expected  that  the  absorption  coefficient  for  type 
II  structures  is  within  an  order  of  magnitude  of  Sample  B. 

Cooley  measured  the  effective  index  of  refraction  of  Sample  B  by  ellipsometry  to 
be  n  =  3.87  (Cooley,  1997:120).  The  closest  approximation  to  the  index  for  the  type  II 
structures  is  based  on  the  index  of  the  GaSb  substrate,  which  is  n  =  4.39  (Madelung, 
1991:120).  The  index  is  required  to  determine  the  reflectivity  of  the  semiconductor,  so 
the  amount  of  power  reflected  can  be  accounted  for.  Without  considering  quantum 
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effects,  39%  of  the  light  traveling  from  air  to  GaSb  will  be  reflected  assuming  normal 
incidence.  The  next  material  in  the  quantum  well  is  AlAsSb,  with  a  lower  index  of 
approximately  n  =  3.2.  Only  0. 1%  of  the  remaining  light  will  be  reflected,  which  is 
negligible.  This  analysis,  however,  does  not  take  into  consideration  the  fact  that  the 
semiconductor  layers  are  all  significantly  less  than  one  wavelength  of  the  incoming  laser. 
Therefore,  experimental  measurements  of  the  reflectance  were  conducted  in  order  to 
reduce  the  uncertainty  introduced  in  Eq  (4.2)  by  the  index  of  refraction.  The  beam  power 
was  measured  before  and  after  reflection  off  each  sample.  The  reflectance  is  then  the 
ratio  of  reflected  to  incident  power.  This  calculation  assumes  a  perfectly  specular 
reflection  at  the  face  of  the  sample.  Eq  (4.3)  was  used  to  determine  effective  indicies  of 
n  =  4.5  for  Sample  B  and  4.3  for  the  type  II  structures  based  on  the  experimental 
reflectance  values. 

Einally,  the  beam  spot  sizes  were  calculated  using  experimental  measurements 
and  Eq  (3.7)  with  an  error  of  less  than  0.2  pm.  The  range  of  spot  sizes  used  in  this 
analysis  accounts  for  the  uncertainty  in  the  calculated  spot  size. 


Table  4.1.  Analysis  of  error  in  unknown  parameters  used  to  calculate  the  carrier  density. 

(1)  determined  experimentally  by  ellipsometry  from  Cooley. 

(2)  determined  experimentally  by  Mode  Master  using  scanning  knife  edge  method. 

Variable 

Value  used 

Range 

%  change  in  carrier 
density  over  range 

a 

2x  10^  cm'^  (1) 

7x  I0^..7x  10^ 

288 

n 

4.5  (type  I) 

4.3  (type  II) 

4.3. ..4.5 

1.8 

2  Wp 

39  pm  (2) 

+/-  0.15  pm 

0.5 

2  Wpi 

58.5  pm  (2) 

+/-  0.20  pm 

<0.1 

39 


Once  the  initial  carrier  densities  are  known  for  various  pump  powers,  a  eurve  ean 
be  fit  to  the  data  to  develop  a  relationship  for  carrier  density  as  a  funetion  of 
luminescenee.  The  initial  earrier  density  is  computed  for  a  number  of  exeitation  powers. 
These  carrier  densities  were  plotted  against  the  experimentally  determined  lumineseenee 
intensities  at  the  same  exeitation  powers.  The  following  fitting  function  was  used  to 
approximate  the  data: 

n(L)  =  uq-L+ui-l'^^  (4  5) 

The  resulting  fit  has  a  R  value  of  0.98  and  is  shown  in  Figure  4.2.  This  fit  is  then 
used  to  ealeulate  all  of  the  carrier  densities  (t  >  0)  from  the  lumineseenee  data  eolleeted. 
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Figure  4.2.  Carrier  density  as  a  function  of  luminescence  intensity.  Data  points  are  the 
initial  (t  =  0)  earrier  densities  generated  by  the  laser  pulse,  which  were  caleulated  using 
Eq(4.4).  The  solid  line  represents  a  fitting  funetion  found  in  Eq  (4.5)  used  to  related 
earrier  densities  at  t  >  0  to  lumineseenee  intensities. 


Curve  fitting  the  carrier  density  rate  equation,  Eq  (4.1),  to  the  data  eolleeted  has 
proven  difficult  (McKay,  2001).  Simultaneously  fitting  four  parameters  that  are  lO"^' 
apart  in  magnitude  is  not  an  easy  task,  and  the  ambiguity  in  reported  recombination 
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coefficients  has  been  noted  (Jang  and  Flatte,  1998).  Previous  researchers  have  used 
different  methods  to  reduce  the  number  of  fitting  parameters  to  solve  this  problem.  The 
work  done  at  the  University  of  Iowa,  cited  in  the  above  paper,  uses  a  calculated  radiative 
coefficient  and  attempts  to  isolate  the  Asrh  coefficient.  McKay  attempted  to  fit  all  three 
coefficients  simultaneously  and  relied  on  the  calculated  initial  carrier  density.  This  paper 
will  use  a  combination  of  these  methods  to  find  the  parameters  that  yield  the  best  fit  to 
the  data. 

SRH  recombination  is  assumed  to  dominate  carrier  recombination  in  the  low 
carrier  regime.  Figure  4.1  shows  carrier  densities  beginning  at  t  =  0  when  the  excitation 
pulse  arrives.  These  points  actually  represent  the  luminescence  decay  13.2  ns  after  the 
previous  excitation  pulse.  The  SRH  coefficient  can  be  estimated  independently  by 
considering  the  carrier  densities  in  the  region  of  2.5  ns  to  13  ns.  If  all  other  terms  in  the 
recombination  rate  equation  Eq  (4.1)  are  ignored,  then  carrier  densities  in  this  region  can 
be  described  by  a  simple  exponential  decay.  A  fit  to  the  data  will  allow  Asrh  to  be 
determined  independently,  reducing  one  of  the  parameters  in  Eq  (4.1).  The  curve  fit  to 
determine  Asrh  for  Sample  91  at  300  mW  excitation  power  with  a  R  =  0.99  is  shown  in 
Eigure  4.3.  The  fit  yields  an  Asrh  =  3.4  x  10  s'  .  SRH  is  commonly  reported  as  a 
lifetime  measured  in  nanoseconds  with  a  value  of  Asrh  The  SRH  lifetime  for  this  fit  is 
29.4  ns. 

The  initial  carrier  density  is  another  key  parameter  required  to  solve  the 
differential  equation  Eq  (4.2).  Although  the  initial  carrier  density  was  calculated 
theoretically,  there  is  error  in  this  term  as  was  discussed  earlier.  If  the  calculated  initial 
carrier  density  is  used  as  the  initial  condition  in  solving  Eq  (4.1),  unwarranted 
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weight  has  been  plaeed  on  this  data  point,  which  forces  the  fitting  curve  to  go  through 
this  point.  The  initial  carrier  density  found  from  Eq  (4.2)  is  used  as  an  estimate  in  the 
iteration  process  used  to  fit  all  four  parameters. 


Figure  4.3.  Curve  fit  to  estimate  Asrh-  SRH  recombination  is  assumed  to  dominate 
after  2.5  ns.  Asrh  is  determined  by  fitting  an  exponential  decay  curve  to  the  data. 

Initial  estimates  for  Asrh  and  the  initial  carrier  density  are  used  in  the 
recombination  rate  equation  Eq  (4.1).  A  range  of  Brad  and  CAuger  coefficients  are  used, 
and  the  resulting  fit  is  compared  to  the  data  set  by  the  following  equation: 

E„or=2(„(q-Yj)^ 

i  ,  (4.6; 

where  n(t)  is  the  predicted  carrier  density  based  on  the  solution  to  Eq  (4.1),  and  Y  is  the 
data  point  taken  at  the  same  time.  A  surface  plot  of  the  squared  error  resulting  from  any 
combination  of  Brad  and  CAuger  coefficients  can  be  made,  as  shown  in  Figure  4.4. 

The  minimum  on  this  surface  represents  the  best  fit  for  a  given  Asrh  and  initial  carrier 
density.  A  number  of  these  surfaces  were  plotted  for  a  series  of  Asrh  and  initial  carrier 
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densities.  While  this  is  a  brute  force  method  of  approaching  the  problem,  it  will  result  in 
the  set  of  coefficients  that  yield  the  smallest  error.  Asrh  and  the  initial  carrier  density 
were  then  varied,  and  the  minimum  error  point  on  the  Brad  and  CAuger  surface  was 
recorded.  The  result  is  another  surface  plot,  shown  in  Figure  4.5,  with  a  minimum  point 
that  represents  the  unique  set  of  coefficients  that  best  fits  the  experimental  data.  This  plot 
was  made  based  on  data  taken  from  Sample  91  at  300  mW.  Each  point  on  the  surface 
represents  the  minimum  error  obtained  from  a  surface  plot  like  the  one  in  Figure  4.4. 


S  Error^ 


Brad 


Figure  4.4.  Error  surface  resulting  from  the  fit  of  Eq  (4.1)  to  the  experimental  data  as  a 
function  of  Brad  and  CAuger-  The  minimum  point  on  this  surface  represents  the  best  fit  for 
a  given  Asrh  and  initial  carrier  density. 
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Asrh 

Figure  4.5.  Error  surface  for  Sample  91  at  300  mW  excitation  power  resulting  from  the 
fit  of  Eq  (4.1)  to  the  experimental  data  as  a  function  of  Asrh  and  initial  carrier  density  no 
Each  point  represents  the  minimum  error  on  the  Brad  and  CAuger  surface  for  a  given  Asrh 
and  no.  The  minimum  point  on  this  surface  represents  the  best  fit  of  the  experimental 
data,  yielding  all  four  fitting  parameters. 


4,1.2,  Recombination  Coefficient  Results 

The  method  describe  above  was  used  to  fit  the  data  taken  from  each  sample.  The 
decay  curve  and  sample  fit  for  Sample  91  are  shown  in  Eigure  4.6.  The  experiment  was 
performed  at  various  excitation  powers.  While  this  causes  the  initial  carrier  density  to 
change,  the  coefficients  that  describe  the  recombination  processes  should  remain  the 
same.  The  Shockley-Read-Hall,  radiative,  and  Auger  coefficients  calculated  for  each 
sample  at  300  mW  are  shown  in  Table  4.2. 
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Table  4.2.  Experimentally  determined  recombination  coefficients  and  initial  carrier 
density  for  each  sample  at  300  mW. 

Sample 

no  (cm'^) 

Asrh'  (ns) 

Brad(cm^  s'^) 

C-Auger  (cm  ) 

B 

7.0+/- 0.1  X  10'“^ 

2  +/-  0.5 

2.7+/- 0.1  X  10'^^ 

3.3+/- 0.1  X  10'^^ 

91 

3.8+/- 0.1  X  lO''^ 

82  +/-  30 

5.8 +/- 0.5  X  10'“ 

6.3  +/-  1  X  10'^^ 

90 

3.5  X  lO*'^ 

60 

2.9  X  10'“ 

5.0  X  10'^^ 

88 

3.6  X  lO*'^ 

39 

3.5  X  10'“ 

6.0  X  10'^^ 

The  error  range  for  Sample  91  resulted  from  the  uncertainty  in  the  absorption  coefficient. 
The  data  analysis  described  in  the  previous  section  was  applied  to  Sample  91  for  both  a 
low  and  high-end  absorption  coefficient  to  study  how  the  uncertainty  in  the  absorption 
affected  the  resulting  coefficients.  The  absorption  coefficient  was  known  for  Sample  B, 
and  the  largest  source  of  error  for  this  sample  is  the  step  size  used  to  vary  each  of  the 
fitting  parameters.  The  error  due  to  background  noise  was  calculated  by  determining  the 
standard  deviation  of  the  luminescence  intensity  when  the  system  sampled  the  same  point 
on  the  decay  curve  for  thirty  iterations.  The  standard  deviation  in  luminescence  intensity 
was  then  converted  to  an  error  in  carrier  density.  It  is  assumed  that  the  background  noise 
will  affect  each  point  the  same.  The  resulting  error  due  to  the  background  noise  was  less 
than  3%  of  the  carrier  density,  and  therefore  was  considered  negligible. 

The  type  II  quantum  well  structures  were  found  to  have  lower  Auger 
recombination  coefficients  than  the  type  I  structure.  This  result  is  in  agreement  with 
previous  research  on  InAs  that  reports  type  II  Auger  coefficients  an  order  of  magnitude 
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lower  than  type  I  Auger  coefficients  (Meyer  et  al,  1998).  The  reduction  in  the  Auger 
recombination  in  this  paper  is  attributed  to  the  removal  of  the  degeneracy  between  the 
band-gap  energy  Eg  and  the  split-off  energy  gap  A  through  band-gap  engineering. 


Carrier 

Density 

(cm'^) 


Figure  4.6.  Carrier  density  decay  over  3-ns  time  period  for  various  excitation  powers. 


4,2,  Carrier  Relaxation 

The  recombination  coefficients  can  be  used  to  quantify  radiative  recombination 
and  nonradiative  recombination  rates,  and  thus  the  efficiency  of  the  device.  Another 
factor  influencing  device  efficiency  is  carrier  relaxation,  or  how  quickly  excited  carriers 
cool  to  the  bottom  of  the  quantum  well.  The  time  resolved  luminescence  data  was  used 
to  study  the  rise  in  intensity  over  the  first  50  ps  after  the  excitation  pulse.  Earlier 
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reporting  indicated  that  most  carriers  have  relaxed  to  the  band  edge  states  within  15-25 
ps  of  excitation  (Jang  et  al,  1996).  The  rise  in  luminescence  intensity  over  time 
corresponds  to  these  carriers  cooling  from  their  initial  energetic  state  high  in  the 
conduction  band  to  the  bottom  of  the  band  where  they  can  recombine  radiatively. 

The  luminescence  rise  was  recorded  for  each  sample  at  various  excitation 
levels,  as  depicted  in  Figure  4.7.  The  rise  time  was  determined  from  the  10%-90% 
intensity  points  on  this  curve.  This  convention  was  previously  used  by  Davis  et  al.  to 
compare  carrier  relaxation  times  in  a  single  quantum  well  structure  (Davis  et  al,  1994). 


Delay  (psec) 

X  300  mW 
•+  225  mW 
1°  150  mW 
'O  75  mW 


Figure  4.7.  Luminescence  rise  time  for  Sample  91  at  various  excitation  powers.  Data 
was  recorded  at  2-ps  intervals.  Zero  on  the  time  scale  represents  the  arrival  of  the 
excitation  pulse. 
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The  subsequent  step  in  their  analysis  was  to  deeonvolve  the  lumineseenee  intensity  data 
with  the  temporal  resolution  of  the  system  in  order  to  aeeurately  measure  the  rise  time. 
Davis  et  al.  used  a  streak  eamera  with  40  ps  resolution  to  obtain  the  time  resolved 
speetra.  This  step  was  omitted  in  the  analysis  here  because  the  temporal  resolution  in  this 
experiment  is  based  on  the  pump  beam  pulse  width  of  122  femtoseconds,  which  will  have 
no  noticeable  affect  on  the  recorded  data.  The  resulting  rise  times  for  the  various  samples 
at  different  excitation  powers  is  shown  in  Table  4.3. 


Table  4.3.  Calculated  rise  time  for  each  sample  at  various  excitation 
powers.  Rise  time  is  measured  from  the  10%-90%  intensity  points. 

Excitation 
Power  (mW) 

Rise  Time  (ps) 

Sample  B 

Sample  91 

Sample  90 

Sample  88 

300 

20 

16 

8 

5 

225 

13 

12 

8 

7 

150 

15 

5 

6 

4 

75 

- 

4 

5 

- 

The  rise  times  for  each  sample  at  varying  excitation  powers  were  also  compared 
graphically.  First,  the  relationship  between  carrier  density  and  luminescence  shown  in 
Section  4.1.1  was  also  used  here  to  covert  luminescence  to  carrier  density.  The  data  was 
normalized  and  fit  to  an  exponential  function  over  the  first  40  ps  to  compare  rise  times,  as 
seen  in  Figure  4.8. 
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Figure  4.8.  Example  fit  with  the  function  shown  to  the  normalized  rise  data  with 
=  0.97. 

The  results  indicated  that  the  rise  time  of  the  luminescence  increases  with  carrier 
density,  as  depicted  in  Figure  4.9  (a).  This  goes  against  a  previous  report  by  Boggess, 
who  found  a  weak  indirect  dependence  of  carrier  density  on  rise  time  (Boggess  et  al., 
2000).  One  explanation  for  the  difference  might  involve  the  density  of  carriers  created  in 
the  sample.  The  initial  carrier  densities  in  this  experiment  are  on  the  order  of  10  cm'  . 
While  Boggess  did  not  report  the  carrier  densities  involved  in  the  cited  work,  previous 
work  by  the  same  group  typically  reported  carrier  densities  on  the  order  of  10  cm' 

(Jang  and  Flatte,  1998;  Flatte  et  al,  1999).  At  higher  carrier  densities  (  n  >  lO'^  cm'^)  , 
hot  phonon  effects  become  observable  (Shah,  1996).  A  non-equilibrium  phonon 
population  arises  as  carrier  cooling  creates  phonons  faster  than  the  phonon  energy  can  be 
dispersed.  The  carriers  can  then  reabsorb  this  phonon  energy,  thus  decreasing  the  carrier 
energy  loss  rate.  Yang  reports  observing  this  effect  in  GaAs  QW  structures  (Yang  et  al., 
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1985).  The  carrier  relaxation  time  is  slower  with  increased  excitation  power,  the  result  of 
which  is  seen  here  in  longer  rise  times. 


Sample  9 1  Rise 


-  225  mW 

—  150  mW 

- 75  mW 

(a)  (b) 


Figure  4.9.  Carrier  Density  rise  time  for  (a)  Sample  91  at  multiple  excitation  powers,  and 
(b)  for  each  sample.  Rise  time  increases  with  carrier  density  and  decreases  with  well 
width  for  the  Type  II  structures. 


Rise  Time  (300mW) 


Delay  (psec) 
Sample  9 1 
Sample  90 
Sample  88 


Figure  4.9  b)  shows  the  rise  times  of  the  different  samples.  The  results  indicate 
that  rise  time  will  decrease  with  the  quantum  well  width  for  the  type  II  structures. 

This  finding  is  in  agreement  with  previous  work  on  the  carrier  relaxation  in 
InGaAs/GaAs  separate  confined  heterostrucutre  lasers.  Davis  et  al.  calculated  the 
multiband  coupling  coefficient,  which  represents  the  wave  function  overlap  integral,  for 
carriers  in  the  quantum  well  and  the  continuum  of  the  sample  tested  (Davis  et  al,  1994). 
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The  coupling  coefficient  is  directly  proportional  to  the  probability  of  a  transition 
occurring,  as  given  by  the  Fermi  Golden  Rule.  The  probability  for  a  transition  between 
states  in  a  quantum  well  was  found  to  be  much  smaller  than  the  probability  of  making  a 
transition  in  bulk  InGaAs/GaAs  (Davis  et  al,  1994).  The  result  is  a  bottleneck  of  carriers 
at  the  highest  energy  state  in  the  quantum  well  as  they  wait  for  the  transition  to  the 
bottom  of  the  well  to  occur.  The  narrower  the  well,  the  less  likely  the  transition  will 
occur,  and  an  increase  in  carrier  relaxation  time  is  observed.  Narrower  quantum  wells 
lasers  will  have  a  higher  threshold  current  density  and  a  lower  overall  efficiency. 

4,3,  Carrier  Temperature 

The  evolution  of  the  carrier  temperature  in  a  semiconductor  is  closely  tied  to  the  rise  in 
the  luminescence  signal  analyzed  in  the  previous  section.  Electrons  initially  excited  by  a 
laser  pulse  will  give  their  energy  to  acoustic  and  optical  phonons,  causing  the  electrons  to 
cool  down  near  the  band  edge  (Shah,  1996:161).  At  high  excitation  powers,  carrier- 
carrier  interactions  dominate,  resulting  in  carrier  temperatures  higher  than  the  lattice 
temperature.  As  time  goes  on,  carrier-optical  phonon  scattering  continues  to  cool 
carriers,  and  the  carrier  temperature  approaches  the  lattice  temperature  (Yang,  1985). 
Auger  recombination  is  dependent  on  the  occupation  of  energy  states  within  the  band 
structure,  which  is  described  by  the  carrier  temperature  (Jang  and  Flatte,  1998). 

4,3,1,  Carrier  Temperature  Results 

Spectral  data  was  recorded  for  Sample  91  and  Sample  B  at  several  times  along  the 
luminescence  decay  curve.  Representative  spectra  for  Sample  91  at  select  time  intervals 
are  shown  in  Figure  4.10.  The  low  energy  side  of  the  curve  is  determined  by  the  band- 
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gap,  as  discussed  in  Chapter  2.  Band  tailing  effeets  eause  lumineseenee  to  be  observed 
below  the  band-gap  energy  in  Figure  4.10.  The  low  energy  side  of  the  eurve  was  fit  to  a 
Gaussian  curve  to  model  line  broadening  and  ensure  it  was  independent  of  time 
(Pankove,  1971:45).  The  empirical  parameter  in  the  Gaussian  fit  has  dimensions  of 
energy  and  deseribes  the  density  of  states,  and  should  remain  eonstant  over  time. 


■  +  500  psee 
O  1000  psee 

Figure  4.10.  Speetral  data  for  Sample  91  at  different  deeay  times.  The  high-energy  tail 
is  fit  with  an  exponential  function  to  determine  earrier  temperature. 


The  Gaussian  fits  to  the  experimental  data  ean  be  seen  in  Figure  4.1 1,  and  the  empirieal 
parameter  is  shown  to  vary  less  than  10%  for  eaeh  eurve  fit. 
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■H  500  psec 
'  O  1000  psec 

Figure  4.1 1.  Gaussian  curve  fits  to  the  low  energy  side  of  the  spectral  data.  The 
empirical  parameter  Eq  listed  describes  the  density  of  states,  and  was  found  to  vary  by 
less  than  10%  over  the  series  of  spectra. 


In  non-equilibrium  situations,  such  as  photoluminescence,  spectral  data  can  be 
used  to  calculate  carrier  temperature.  The  high-energy  portion  of  the  luminescence  curve 
is  dictated  by  the  Fermi  distribution  function  described  in  Chapter  2. 

If  the  carriers  have  a  large  amount  of  energy,  as  is  the  case  in  this  experiment,  the  Fermi 
function  can  be  approximated  by: 

-E 

f(E)  =  Ae^'^  ^  (4  7) 
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which  is  the  Bolzmann  distribution  function,  where  T  is  the  carrier  temperature.  Eq  (4.7) 
can  be  used  to  fit  the  high-energy  side  of  the  luminescence  curve,  and  the  carrier 
temperature  at  the  time  the  speetra  was  taken  can  be  determined,  as  shown  in  Figure  4.12. 


Energy  (eV) 

50  psec 
+  550  psec 
^  950  psec 

^  1550  nsec 

Figure  4.12.  Boltzamn  distribution  curve’s  used  to  extract  carrier  temperature  shown  in 
figure.  Average  carrier  temperature’s  cool  as  relax  to  the  band  edge  and  recombine. 
This  method  is  applied  to  spectra  taken  at  different  times  along  the  lumineseenee  decay 
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curve  to  determine  how  earrier  temperature  evolves  over  time.  This  method  was  first 
used  in  non-equilibrium  situations  such  as  this  by  Linde  to  measure  the  change  in  carrier 
temperature  over  time  in  bulk  GaAs  (Linde,  1978).  Jang  et  al.  have  reported  carrier 
temperature  evolution  in  mid-IR  MQW  structures  using  exponential  fitting  functions 
(Jang  eta/.,  1997). 


4,3,2.  Carrier  Temperature  Modeling 

The  time  dependence  of  carrier  temperature  can  be  modeled  by  examining  the 
carrier  energy-loss  rate  with  a  method  described  by  Shah  (Shah,  1996:162-165).  The 
average  energy-loss  rate  per  carrier  to  the  lattice  can  be  described  by: 


dt 


-E 


^  3D 


(4.8) 


where  E  is  the  longitudinal  optieal  (LO)  phonon  energy,  and  T  is  the  carrier  temperature 
(Shah,  1996:163).  Optical  phonon  scattering  is  dominant  at  carrier  temperatures  greater 
than  40  K,  and  will  be  considered  the  only  source  of  energy  loss  in  this  analysis  (Y ang, 
1985).  l/x3D  is  the  characteristic  rate  for  polar  optical  phonon  scattering  and  is 
determined  by  (Shah,  1996:163): 


-1 


2 

e  -71 


2-mE 


h 


^  sj 


(4.9) 


where  m  is  the  carrier  mass,  and  Si  and  Ss  are  optical  and  static  dielectric  constants.  The 
carrier  mass  is  estimated  to  be  m  =  2.095  x  10'  kg,  the  effective  electron  mass  in  InAs  at 
300  K  (Madelung,  1991).  The  optical  and  static  dielectric  constants  for  InAs  are  12.25  So 
and  15.15  So,  respectively  (Madelung,  1991). 
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The  optical  phonon  scattering  rate  in  Eq  (4.9)  applies  to  bulk  semiconductors  and  is 
related  to  quantum  well  structures  by: 

2 

2D  =  ~  3D 

^  .  (4.10) 

The  optical  phonon  energy  is  taken  from  literature  to  be  approximately  28  meV 
for  mid-IR  semiconductors  (Jang  et  al.,  1996;  Marciniak,  1995).  Using  these  values  for 
an  InAs  quantum  well  structure,  the  longitudinal  phonon  scattering  lifetime  is  estimated 
to  be  143  fs.  The  average  energy-loss  rate  by  a  carrier  to  the  lattice  can  now  be  described 
as  a  function  of  temperature,  as  seen  in  Figure  4.13. 


Figure  4.13.  Calculated  optical  phonon  energy-loss  rate  as  a  function  of  carrier 
temperature  for  InAs.  At  low  temperatures,  acoustic  phonons  dominate  resulting 
in  higher  loss  rates  below  40  K,  which  is  not  depicted  in  this  figure. 
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This  average  energy  loss  rate  ean  then  be  equated  to  the  average  loss  rate  of 
temperature  by  the  following  equations  (Shah,  1996:165): 

^  d- 
dt  dt 

9 

where  the  average  energy  E  =  3kTc,  and  Tc,  and  To  are  the  earrier  and  initial 
temperatures.  Eq  (4.1 1)  is  numerieally  inverted  to  deseribe  the  earrier  temperature  deeay 
as  a  funetion  of  time. 

While  the  model  used  to  determine  energy-loss  rates  in  Eq  (4.1 1)  is  useful  in 
determining  how  energy  is  transferred  from  earriers  to  phonons,  several  assumptions 
were  made  that  will  eause  experimental  results  to  differ  from  theory.  Two  notable 
assumptions  are:  1)  phonon  oeeupation  remains  unaffeeted  by  the  photoexeitation,  and  2) 
sereening  affeets  have  been  ignored  (Shah,  1996:171-173).  The  danger  in  making  the 
first  assumption  is  that  earners  ean  reabsorb  phonons  if  the  number  of  phonons  is  greatly 
inereased  by  the  exeitation  pulse,  eausing  less  energy  to  be  lost  than  predieted.  This  is 
the  hot  phonon  effeet  diseussed  in  Seetion  4.2.  The  seeond  assumption  ignores  the 
eomplex  interaetions  between  earriers. 

As  diseussed  in  Seetion  2.2,  the  peak  lumineseenee  from  an  exeited  sample  shifts 
to  lower  energies  as  eleetrons  eool  from  their  exeited  state.  The  high-energy  side  of  the 
speetrum  is  fit  with  the  exponential  in  Eq  (4.7)  to  determine  earrier  temperature  as 
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discussed  in  the  previous  section.  The  temperature  results  are  now  plotted  in  Figure  4.14 
along  with  the  theoretical  expression  for  temperature  derived  in  Eq  (4.1 1).  A  fitting 
parameter  C  =  0.5  is  used  to  compare  the  model  carrier  temperature  to  experimental 
results. 


Carrier  Cooling  Curve 


-  Theoiy  Time  (psec) 

X  Sample  91 
O  Sample  B 


Figure  4.14.  Carrier  temperature  as  a  funetion  of  time  for  Sample  91  and  Sample  B.  The 
theoretieal  model  for  carrier  temperatures  in  InAs  is  plotted  as  a  solid  line.  An  initial 
temperature  of  350  K  and  a  fitting  parameter  C  =  0.5  was  used  to  fit  the  data. 


This  fitting  parameter  indicates  that  the  carrier  loss  rate  observed  is  half  the  rate 
predicted  by  the  model.  Factors  that  reduce  carrier  loss  rates  include  the  reabsorption  of 
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phonons  and  carrier  screening  discussed  in  the  previous  section.  Fitting  parameters 
between  0.1  -  0.5  have  been  used  to  fit  previous  experimental  data  to  this  model  for 
GaAs  and  InGaAs  (Leheny  et  al., 1919;  Kash  et  al,  1984). 

The  carrier  temperature  data  in  Figure  4.13  indicates  that  carriers  initially  cool 
very  rapidly  within  the  first  400  ps,  and  then  the  cooling  rate  begins  to  taper  off.  At  this 
point  the  average  carrier  energy  has  fallen  far  below  the  characteristic  optical  phonon 
energy,  which  is  the  prime  source  of  carrier  cooling  (Jang  et  al,  1996).  The  carrier 
temperatures  have  not  reached  the  lattice  temperature  (T  =  77  K)  in  the  time  scale 
observed.  This  is  likely  a  result  of  Auger-induced  heating,  which  contributes  energy 
comparable  to  the  band-gap  energy  every  time  an  Auger  recombination  occurs  (Jang  and 
Flatte,  1998).  Carrier  temperatures  in  Sample  B  appear  to  be  at  least  initially  higher  than 
Sample  91,  which  also  has  a  lower  Auger  coefficient.  This  reiterates  the  point  that  carrier 
cooling  and  recombination  processes  are  interrelated  events  that  dictate  the  efficiency  of 
a  semiconductor  device. 


59 


5,  Conclusion 


The  objective  of  this  research  was  to  study  how  different  semiconductor 
structures  impact  carrier  recombination  and  carrier  relaxation.  The  goal  is  to  provide 
crystal  growers  feedback  on  device  structures  in  an  effort  to  increase  the  overall 
efficiency  of  mid-IR  semiconductor  lasers.  This  chapter  will  summarize  key  findings  and 
make  suggestions  for  future  work. 

5.1,  Summary 

Data  from  the  TRPL  experiment  was  used  to  investigate  the  recombination 
coefficients,  carrier  relaxation,  and  carrier  cooling.  The  recombination  coefficients  were 
calculated  by  fitting  the  recombination  rate  equation  to  carrier  densities  over  time.  Two 
significant  sources  of  error  in  determining  the  coefficients  are  highlighted.  First,  the 
absorption  coefficient  of  the  semiconductor  leads  to  the  greatest  error  in  computing  the 
initial  carrier  density.  Second,  the  fitting  process  used  in  previous  work  to  calculate  the 
coefficients  did  not  vary  every  parameter  in  the  recombination  rate  equation.  To  find  the 
parameters  that  best  fit  the  experimental  data,  each  coefficient  and  the  initial  condition 
was  varied  to  minimize  the  error  in  the  fit.  The  type  I  structure  was  found  to  have 
significantly  higher  Auger  coefficients  compared  to  the  type  II  structures. 

Analysis  of  the  carrier  rise  time  indicates  that  the  carrier  relaxation  time  will 
increase  at  higher  excitation  powers  and  with  narrower  wells.  Hot  carrier  effects  cause 
carriers  to  reabsorb  phonon  energy  and  increase  their  relaxation  time.  Narrower  wells 
have  a  lower  probably  of  downward  transitions  to  the  band  edge  than  wider  wells,  which 
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causes  an  increase  in  relaxation  time  due  to  the  carrier  bottleneck  that  occurs  at  the  top  of 
the  narrow  well. 

Spectral  data  from  the  TRPL  experiment  was  used  to  extract  the  carrier 
temperature  over  time.  It  was  found  that  type  I  structures  had  higher  carrier  temperatures 
than  type  II  structures,  which  corresponds  to  the  Auger  coefficients  found  for  each 
structure.  The  Auger  process  adds  energy  to  the  system  by  promoting  a  carrier,  which  is 
measured  as  warmer  carrier  temperatures.  These  results  confirm  that  type  II  structures 
have  reduced  the  impact  Auger  recombination  in  mid-IR  semiconductor  devices. 

5,2,  Future  Work 

Future  work  on  this  project  should  be  focused  in  two  areas.  First,  improvements 
can  be  made  in  the  analysis  to  increase  the  accuracy  of  the  recombination  coefficients. 
Additionally,  the  TRPL  experiment  can  be  converted  with  some  effort  to  examine  mid-IR 
semiconductor  lasers. 

As  cited  earlier.  Auger  recombination  coefficients  reported  here  and  elsewhere 
vary  over  several  orders  of  magnitude.  To  accurately  compare  the  efficiency  of  different 
structures,  the  source  of  these  fluctuations  must  be  identified.  This  research  indicates  that 
the  absorption  coefficient  plays  a  significant  role  in  determining  the  initial  carrier 
densities.  Absorption  coefficients  should  be  measured  for  each  sample  through 
ellipsometry  or  other  means  to  reduce  the  error  in  the  initial  carrier  density.  In  addition, 
the  curve  fitting  process  described  in  this  work  could  be  further  automated,  which  would 
provide  more  timely  results  and  allow  for  further  analysis. 
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One  important  note  in  the  eolleetion  of  speetral  data  is  that  the  upeonversion 
efficieney  remained  high  (  >  90%)  when  performing  the  experiment  over  the  speetra  of 
Sample  91.  This  indieates  that  it  may  not  be  neeessary  to  eonstantly  rotate  the  non-linear 
erystal  in  order  to  reeord  speetral  data  as  suggested  in  earlier  work  (Franz,  1997). 

The  method  used  in  this  experiment  to  eolleet  time-resolved  speetra  has  broader 
applieations  than  determining  earrier  temperatures.  Laser  engineers  at  Lineoln 
Laboratory  at  the  Massachusetts  Institute  of  Technology  have  theorized  that  unintended 
modes  may  exist  on  a  very  short  time  scale  in  the  output  of  semiconductor  lasers.  These 
modes  will  reduce  the  gain,  and  therefore  the  output  power  in  the  desired  mode.  The 
TRPL  experiment  could  be  converted  to  examine  the  output  from  an  optically  pumped, 
mid-IR  semiconductor  laser  to  test  this  theory.  The  sample  laser  devices  provided  by 
MIT/LL  require  that  the  pump  beam  provide  optical  confinement  in  the  structure.  This 
requires  that  the  pump  beam  must  be  focused  on  to  the  sample  with  a  cylindrical  lens 
resulting  in  a  focused  spot  on  the  order  of  1  mm  by  100  microns.  The  pump  beam  must 
also  excite  the  sample  from  the  side  rather  than  from  the  front  as  depicted  in  Figure  3.3. 
With  these  exceptions,  the  TRPL  experiment  is  well  suited  to  investigate  laser  devices. 
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Appendix  A. 


A.1.  Crystal  Rotation  Angle 


Appendix  A.  1.  is  a  program 
derived  from  Lt  Col  Marciniak's 
Mathmatica  program  deigned  to 
calculate  crystal  rotation  based 
on  PL  wavelength. 


KTA  Rotation  Angle 


Index  ellipsoid  for  KTA  crystal 


iix(x) 


1  23522  2 

1.90713  + - — - .01025  X 


1  _  (^->9692^ 

X  J 


n^x):= 


2.15912  + - L22222 - .01096- 


1  - 


.21844^ 

>.  J 


•- 


2.14786  + - L??559 - .oi436-A.^ 


1  - 


.22719^ 

X  J 


X:=.3,.3l..4  A.  in  pm 
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Inputs:  X.pL  ;=  3.5  A.  in  fim 


^Pump  •- 


^SFG  := 


1 

1  ^  1 
^PL  ^Pump 


X-SFG  =  0.657113 


n^Xph)  =  1 .740622  nx(A,Pump)  =  1 .792606  n^Xpump)  =  1 .8828 1 1 


ny(^SFG)=  1-811004 


Determining  extraordinary  index 


ne(0Pump)  ' 


%(^Pump)'*^^Pump) 


^nx( A-Pump)  ‘  cos|^0  pump‘  ^  +  ^*'z{ A-Pump) '  sm^O  Pump'  j 


ne(46.49324)  =  1.833707 
6  :=0..90 


Snell's  law 


Crystal  cut  at  43  degrees 
20  degrees  between  6pl  and  Opump 


e 
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180 


sin 

n 

180 

180  • 

20 - asm 

n 

**e(0Pump)'®™ 

(47-epump)-  jgQ 

n>(^PL) 

+  47 


0SFGl(0Pump)  :=  —  asin 

71 

nj(A.pL)sin|^0PL(0Pump)~j 

ae(0  Pump)  Pump' 

^SFG 

^  ^PL 

+ 

^Pump  j 

ny(^SFG)_ 

0SFG2(0Pump)  :=  — acos 
n 

Y  (  7C  ^ 

nj(>.PL)  cos[  0PL(0Pump)  — 

ne(0  Pump)  ■ ‘XJS^0  Pump' 

^SFG 

^  ^PL 

•  + 

^Pump 

)  ny(^SFG)_ 

Guess  value  for  intersection 

9  Pump  •=  40 

Given 

eSFGl(0Pump)  =  0SFG2(9pump) 
X  :=  Find(epump) 

x=  46.493242 
0Pump  •=  X 

0PL(x)  =  57.819271 

0SFGl(x)  =  48.524166 
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6R0T(ePump)  :=  -^•asinj^ny(A,pL)-sinj^(epL(epump)  - 


0ROT(x)  =  19.070727 


0SFGout(0Pump)  :=  —  asin  nj(A.SFG)  sin  (0SFGl(0Pump)  -  47)--^ 

7C  1  oO 


0SFGout(x)  =  2.761013 


A.2.  Finding  the  Phase  Mismatch 


Calculate  Ak 

Phased  matched  condition 


Akx:= 


ny(^PL)cos|^0PL(0pump)-^j  ne(0pump)cos|^0Punip-^j  «v(^SFG)-cos^0SFGl(0Pump)-^ 


^PL 


^Pump 


^SFG 


Ak, 


•^(^Pl)  .  \  X  ^  ne(0Pump)  .(  jt  ^  nj(^SFG)  .[  ,  ^  x 


+  Akz 


Ak=  8.836879  X  10 


Close  to  zero,  as 
expected 


The  PL  wavelength  will  be 
changed  to  determine  Ak 


X.pL  :=  3.8 


0Pump:=  46.49324 


^SFG  := 


1 

1  1 
^PL  ^Pump 


XsFG  =  0.666999 


180  . 

sin 

n 

180 

180  . 

20 - asm 

n 

ne(0Pump)sin 

'  ■  ■■“  *3sin 

7t 

ny(^PL) 

180 


+  47 
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180 

n 


ny(XpL)-sin|^0PL(epump)-^^ 

1  Pe(0Pump)-®'*’^0Pump' 

^SFG 

1  ^PL 

"  + 

^Pump  y 

ny(^SFG)_ 

ny(XpL)cos|^epL(epump)— j  ne(0  Pump) -cosine  Pump- ny(A,SFG)-cos^0SFGl(ePump)-^ 

^PL  ^Pump  ^SFG 


Ak, 


")(^Pl)  .  \  "e(0Pump)  .  7t  •^^SFg)  .  ^  n 


Ak 


.  2 

:-  ^  Akx  + 


Aky 


Ak=  2.968829  X  10 


-3 


Need  to  determine  effective  length 

Quasi-statci  interaction  length 


t:=  122-10 

ne  :=  1 .833707  no  :=  1 .740662 
c:=310* 


-  10 


o  =  3.1015  X  10 


AX  :=  10-10  ’  X:=  809-10  ’  A©  := 


2-7t-c-AX 


Lqs  :=  (u-A©)  Lqs  =  1.1 195  x  10 
L:=112  urn 


-4 


The  efficiency  equation  can  now  be  used 

^2 


T1  := 


sin(L-Ak) 


(L-Ak) 


Ti  =  0.963685 
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The  upconversion  efficiency  is  96.4%  of  the  peak  efficiency  if  upconverting  3.8  pm  PL  with  the 
crystal  tilt  optimized  for  3.5pm  PL. 


Efficiency  Curve  for  A.up  near 
^peak  =  3.5pm 


A.3.  Carrier  Temperature 

Constants 

h:=  6.6260755- 10“^'* 
k;=  1.38065-10“^^ 

q;=  1.602177-10“*’ 

E  :=  28-10“  ^-q  Energy  of  optical  phonon  taken  from  Jang  1 996 

Eo:=  8.854187817-10“ 
ei  :=  12.25-eo 

These  are  dielectric  permittivities  taken  for  InAs 

Es  15.15-Eo 

m  :=  .023-9.1093897-10”  Electron  mass  for  InAs 

1 

^3D  :=  — — ; - - 
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t  3l)  =  2,250005  X  10 


13 


t2n  ;=  — ^30 

-E 

H  k-T 
ilF<T);= - c 

''2Dq 
T  :=  1  „  KKK) 


X2D  =  1.432398  X  10 


-  13 


10 


dl;(  1 50)  =  2.240455  x  10  22  meV  /  ps 


-.vin 


l(l'c)  :=  (I'c  -  To)'  '  c  To 0.  1000 


C<xiling  Ciir\'e 


Sample  cooling  curve  for 
electrons  in  InAs  with 
initial  temp  =  750  and 
C  =  0.5 


69 


Appendix  B 


B,l.  Decay  Curves 


Carrier  Density 
(em'^) 


Sample  B 


Time  tsect 
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3 

Carrier  Density 
(em'^) 

2 

1 


Carrier  Density 
(em'^) 


Sample  88 


Time  (see) 


Intensity  (a.u.) 


B,2,  Raise  Time  Experimental  Data 


2.5-10 


2-10 

^  4 

^  1.5-10 

^  A 

4 

I  1-10 


Sample  91  Rise 


5000 


T 

>< 


T 


xx^x^x  ^^x 

jdin^n^nnnnn^nni^i 

X  oO°ooo°OoO°oOooO( 


■10 


10 


20 


30 


40 


Delay  (psec) 

X  300  mW 
+  225  mW 
□  150  mW 
O  75  mW 


Delay  (psec) 
X  300  mW 
+  225  mW 
□  150  mW 
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Intensity  (a.u.) 


Delay  (psec) 
X  300  mW 
-H  225  mW 
n  150  mW 
O  75  mW 


-I-  225  mW 
O  150  mW 
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Intensity  (a.u.)  Intensity  (a.u.) 


B,3,  Rise  Time  Curve  Fits 


Sample  88  Rise 


225  mW 
150  mW 


Sample  91  Rise 


Delay  (psec) 
300  mW 
225  mW 
150  mW 
75  mW 


Sample  90  Rise 


225  mW 
150  mW 


Sample  B  Rise 


Delay  (psec) 
300  mW 
225  mW 
150  mW 
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without  emitting  a  photon.  Nonradiative  recombination  mechanisms  occur  over  a  short  time  period  and  difficult  to  measure. 
Growth  methods  have  significantly  reduced  the  nonradiative  recombination  in  some  materials.  The  objective  of  this  research 
is  to  further  the  understanding  of  how  quantum  well  structures  impact  carrier  recombination.  InAsSb/InAlASb  and 
InAs/GalnSb  quantum  well  structures  were  studied  with  time-resolved  photoluminescence  utilizing  upconversion,  a  non-linear 
wave  mixing  technique.  This  research  reports  Shockley-Read-Hall,  radiative,  and  Auger  recombination  coefficients  at  77  K. 
The  luminescence  rise  times  of  type  I  and  type  II  structures  are  also  compared.  The  number  of  states  available  within  the 
quantum  well  was  found  to  dictate  how  quickly  carriers  were  able  to  recombine  radiatively.  Finally,  spectral  data  was  taken 
to  examine  the  spectral  decay  of  the  luminescence.  Carrier  temperatures  were  extracted  from  the  spectral  data.  Type  I 
structures  were  found  to  have  hotter  carrier  temperatures  and  higher  Auger  coefficients  than  type  II  structures. 
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